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1.0  Introduction 


2 

Integrated  injection  logic  (I  L)  gates  are  a class  of 

integrated  circuits  which  can  be  fabricated  on  bipolar 

transistor  production  lines,  making  it  possible  to  combine 

2 

other  bipolar  circuitry  with  I L logic  on  the  same  inte- 
grated circuit  (IC)  chip.  They  promise  high  packing  den- 
sities and  low  delay-power  products,  making  them  potentially 
competitive  with  conventional  metal-oxide-semiconductor 
(MOS)  logic  circuitry.  The  fundamental  I L logic  circuit 
is  a simple  inverter  (see  sketch)  which  physically  consists 
of  a vertical  npn  multiple  emitter  transistor  operated  in 
the  inverse  mode  so  that  the  conventional  emitters  perform 
as  collectors.  Base  drive  to  the  npn  inverter  is  supplied 
by  a lateral  pnp  transistor  whose  p-type  collector  is  inte- 
grated (or  merged)  with  p-type  base  of  the  inverter. 


The  intended  end  use  of  this  mathematical  model  places 

limitations  upon  its  complexity.  Physical  mechanisms  of 
2 

I L operation  must  be  approximated  with  the  degree  of  ac- 
curacy needed  for  engineering  purposes,  yet  the  model  must 

remain  simple  from  a computational  point  of  view.  During 
2 

I L design  the  engineer  should  be  able  to  undertake  an 
extensive  program  of  computational  experimentation,  with- 
out involving  an  excessive  amount  of  computer  time. 

The  initial  plan  for  this  project  was  to  complete 
the  development  effort  over  a two-year  period.  The  first 
year  was  to  be  directed  toward  modeling  the  vertical  tran- 
sistor  in  an  I L structure,  without  including  the  mechanisms 
of  radiation  damage.  Thereafter,  during  the  second  year, 
it  was  planned  to  complete  this  model  development  effort. 

During  the  first  year  it  was  also  planned  to  develop 
a two-dimensional  mathematical  model  for  vertical  transis- 
tor operation.  For  practical  reasons,  this  two-dimensional 

2 

model  was  not  intended  as  an  I L design  tool  but,  instead, 

as  a research  tool  to  be  used  during  the  development  of  a 

simplified  engineering  model.  Clearly,  the  funding  level 

for  this  program  would  not  permit  the  complete  development 

of  such  a two-dimensional  computer  model.  Therefore,  it 

was  our  intention  to  modify  an  existing  model  for  bipolar 

transistor  operation;  modifications  that  would  permit  an 

2 

analysis  of  device  operation,  as  encountered  in  the  I L 
vertical  transistor. 
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It  can  be  stated  that  these  modifications  far  exceeded 


k 

s 


i 


our  expectations.  At  the  initiation  of  this  program  a two- 
dimensional  model  for  the  bipolar  transistor  was  opera- 
tional; it  was  not  a good  computer  model,  but  it  worked. 

After  introducing  the  required  modifications  many  difficul- 
ties were  encountered.  These  difficulties  were  primarily 
associated  with  computational  stability.  Algorithms  that 

worked  for  conventional  bipolar  transistor  operation  proved 
2 

inadequate  for  I L operation,  and  many  of  these  difficulties 

could  not  be  overcome.  Nonetheless,  this  difficulty  did 

not  alter  the  program  of  development  directed  toward  an 

engineering  model  for  this  semiconductor  structure. 

2 

The  engineering  model  proposed  here  for  I L operation 
is  not  new.  In  fact,  this  model  represents  a new  application 
of  models  previously  developed  for  bipolar  transistor  opera- 
tion [1,2],  As  a consequence,  a large  part  of  this  pro- 
posed model  has  been  used  extensively  for  device  design, 
and  it  has  proven  adequate  for  engineering  purposes.  Fur- 
ther, it  is  known  that  this  model  is  computationally  fast; 

a suitable  implementation  of  this  model  would  permit  the 
2 

analysis  of  an  I L vertical  transistor  in  an  estimated  20 

to  30  seconds  of  CPU  time  (assuming  an  IBM  360/85) . 

Described  here  is  an  implementation  of  this  model,  and 

2 

Its  application  to  specific  problems  of  I L device  design. 


[1]  D.P.  Kennedy  and  P.C.  Murley,  IBM  Jour.  , 8^,  482  (1964). 

[2]  D.P.  Kennedy  and  P.C.  Murley,  Solid-State  Electronics, 
15,  203  (1972). 
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The  discussions  presented  in  this  report  are  limited  to 
those  aspects  of  model  development  completed  during  the 
first  contract  year.  Because  this  program  of  effort  is 
not  being  funded  to  its  completion,  further  developments 
of  this  model  will  be  described  elsewhere. 

2 

2.0  Mathematical  Model  for  I L Operation 

Figure  1 illustrates  the  vertical  transistor  assumed 

in  this  proposed  model  for  the  vertical  transistor  of  an 
2 

I L type  structure.  Simplification  of  this  model  Is 

accomplished  by  its  subdivision  into  two  fundamer  ally 

different  regions;  an  intrinsic  region  and  an  extrinsic 

region.  Fig.  2.  Clearly,  this  simplification  ignores  two 

2 

important  aspects  of  I L operation.  First,  this  simpli- 
fication ignores  mechanisms  encountered  at  the  up  emitter 
periphery,  and,  second,  it  ignores  the  two-dimensional 
flow  of  electric  current  within  the  transistor  base  region. 

Peripheral  mechanisms  associated  with  this  structure 
are  recognized  as  important  to  vertical  transistor  opera- 
tion. Further,  it  is  recognized  that  these  mechanisms 
are  greatly  influenced  by  the  close  proximity  of  the  injec- 
tor junction.  For  this  reason,  it  was  [lanned  that  the 
model  for  these  peripheral  mechanisms  would  be  included  in 
the  second  year  effort,  as  a part  of  the  horizontal  transis- 
tor. Despite  the  absence  of  these  mechanisms,  it  is  shown 

2 

that  in  its  present  form  this  model  for  a vertical  I L 
transistor  provides  extensive  insight  and  quantitative 
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information  concerning  the  operation  of  such  a semi- 
conductor device. 

During  up  operation,  base  region  electric  current 
produces  a two-dimensional  potential  distribution  within 
this  semiconductor  structure.  This  potential  distribution 
is  known  to  produce  a variation  of  forward  bias  across 
various  regions  of  the  up  emitter  junction.  Thereby,  the 
apparent  current  gain  under  up  operation  will  differ  be- 
tween each  collector  junction;  this  gain  will  decrease  with 
an  increase  of  distance  from  the  base  region  electrical 
contact  and/or  the  injector  junction. 

It  is  suggested  that  this  situation  is  similar,  in 

many  ways,  to  conventional  bipolar  transistor  operation 

when  driven  into  saturation.  Further,  it  is  suggested  that 

2 

this  mechanism  can  be  adequately  modeled  for  I L operation 
using  techniques  similar  to  those  used  for  the  bipolar 
transistor.  Specifically,  lumped  model  approximations  for 
each  vertical  transistor  region  could  be  coupled  together  by 
resistors  approximating  the  base  resistance  of  this  struc- 
ture. Although  this  technique  would  complicate  the  present 
model,  it  could  provide  a means  to  approximate  variations  of 
current  gain  due  to  the  two-dimensional  base  region  poten- 
tial distribution. 

2 . 1 The  Physical  Model 

The  physical  model  for  this  semiconductor  device  must 
be  representable  in  terms  of  a specified  impurity  atom 


profile.  Therein  we  have  a topic  of  much  discussion 
throughout  the  semiconductor  industry.  It  is  frequently 
assumed  that  an  inability  to  adequately  approximate  these 
impurity  atom  distributions  limits  the  validity  of  models 
for  bipolar  transistor  operation.  It  is  proposed  that 
this  inability  produces  only  a need  to  understand  these 
inherent  differences  between  experiment  and  theory.  Such 
understanding,  once  attained,  permits  modification  of  the 
mathematical  model  in  a fashion  consistent  with  experiment. 

It  is  recognized  that  few  semiconductor  devices  contain 
impurity  atom  distributions  that  are  ideally  represented  by 
closed  mathematical  expressions  of  the  type  derived  in  text- 
books on  this  subject.  Nonetheless,  it  is  suggested  that 
these  experimental  anomalies  appear  primarily  in  regions 
containing  an  exceedingly  large  impurity  atom  density.  It 
is  also  suggested  that  tradiational  closed  mathematical  ex- 
pressions can  adequately  approximate  the  active  regions  of 
a bipolar  transistor — assuming  these  expressions  are  modified 
to  correct  for  the  high  concentration  anomalies.  Here  we 
have  the  basis  for  our  model  of  the  impurity  atom  distribu- 
tions within  the  intrinsic  and  extrinsic  regions  of  a bipolar 
transistor. 

2 

Specifically,  it  is  assumed  the  base  region  of  an  I L 
vertical  transistor  is  p-type,  and  it  has  been  fabricated 
by  diffusing  boron  into  silicon.  Further,  it  is  assumed 
that  this  diffusion  process  is  of  the  traditional  two-step 
type — predeposition  and  oxidation.  It  has  been  shown  [3] 

[3]  D.P.  Kennedy  and  P.C.  Hurley,  Proc.  IEEE,  52,  372  (1963). 
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that  the  oxidation  phase  of  this  diffusion  process  yields 
an  impurity  atom  distribution  that  is  well  approximated 
by  a Gaussian  type  of  distribution, 


C(x) 


Co6xp 


Ck+X] 

2 

i 2L  J 

(1) 


where  represents  the  impurity  atom  surface  concentration, 

X is  the  distance  from  the  silicon-oxide  interface,  K is 
the  oxide  thickness,  and  L is  the  impurity  atom  diffusion 
length. 

In  the  derivation  of  Eq.  1 we  neglect  changes  of  im- 
purity atom  distribution  arising  from  the  dopant  solubility 
in  Si02.  Because  this  mechanism  introduces  a relatively 
short-range  modification  of  the  impurity  atom  distribution, 
it  has  been  found  to  have  small  influence  upon  bipolar 
transistor  operation.  For  this  reason,  mechanisms  influencing 
the  impurity  atom  distribution  near  the  oxide-semiconductor 
interface  are  not  considered  in  this  analysis. 

Due  to  the  initial  predeposition  phase  (high  impurity 
atom  concentration)  it  has  been  found  that  the  impurity 
atom  diffusion  length,  L,  is  seldom  known,  a priori.  Each 
different  diffusion  process  produces  a different  diffusion 
length  and,  as  a consequence,  we  seldom  know  C (x)  at  any 
given  location.  Despite  this  situation,  when  boron  is 
diffused  into  n-type  silicon  of  known  background  doping, 

Cg,  we  have 
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For  the  analysis  of  semiconductor  devices  we  need  only  know 
where  the  junction  is  located,  and  the  surface  concentration 
of  boron,  C^,  which  can  be  determined  experimentally. 

It  has  been  found  that  Eq.  2 adequately  approximates 
the  boron  impurity  atom  distribution  throughout  the  active 
regions  of  a bipolar  transistor.  This  verification  is  in- 
direct, but  extensive.  Specifically,  numerous  device  para- 
meters are  exceedingly  dependent  upon  the  base  region  im- 
purity atom  distribution;  for  example,  the  reach- through 
voltages.  Satisfactory  agreement  is  attained  with  experiment 
when  these  parameter  calculations  are  based  upon  a base  region 
impurity  profile  of  the  type  given  by  Eq.  2. 

Approximating  the  down  emitter  impurity  profile  in 
this  device  represents  a difficult  problem  of  judgment. 
Specifically,  the  narrow  base  width  used  in  today's  technology 
makes  it  necessary  to  use  a high  concentration  arsenic  dopant 
to  form  the  down  emitter.  This  becomes  necessary  to  prevent 
compensation  of  the  base  region  impurities  and,  thereby, 
unduly  reduce  the  overall  base  region  impurity  atom  density. 
Phosphorous,  for  example,  produces  a substantially  smaller 
impurity  atom  gradient  than  does  arsenic,  and  phosphorous 
would  restrict  I L structures  to  a relatively  large  base  width 
Contrasting  with  phosphorous,  the  arsenic  diffusing  process 
produces  an  impurity  atom  profile  for  which  there  is  no 


available  mathematical  model  [4];  the  resulting  profile 
is  exceedingly  steep  in  the  vicinity  of  the  collector 
junction,  Fig.  3.  As  a consequence,  any  mathematical  model 
of  this  diffusion  process  must  be  an  approximation  that  is 
consistent  with  experiment,  yet  not  alter  the  electrical 
properties  of  any  device  under  consideration. 

Experience  has  shown  that  this  arsenic  profile  can  be 
adequately  approximated  by  a complementary  error  function, 
using  an  unrealistically  large  impurity  atom  surface  con- 
centration. For  example,  throughout  this  calculation  we 

24 

use  an  assumed  surface  concentration  of  about  10  atoms/ 

3 24 

cm  . This  is  not  intended  to  imply  that  we  believe  10 

3 

atoms/cm  represents  a realistic  value;  it  is  obviously 
unrealistically  large.  Nevertheless,  by  this  means  we 
attain  the  needed  impurity  atom  gradient  near  the  up  collector 
junction  and,  thereby,  produce  a realistic  amount  of  base 
region  impurity  atom  compensation.  This  unrealistic  value 
of  surface  concentration  yields  near  the  down  emitter  junc- 
tion an  impurity  atom  gradient  typical  of  that  produced  by 
arsenic.  Further,  this  approximation  technique  yields  a 
similar  type  of  impurity  distribution  throughout  regions 
adjacent  to  this  junction,  regions  that  are  associated  with 
transistor  operation. 

A similar  approach  has  been  taken  in  our  approximation 
of  the  buried  layer.  Specifically,  throughout  this  analysis 

[4]  D.P.  Kennedy  and  P.C.  Murley,  Proc.  of  IEEE,  59,  335 

(1971)  . 
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we  assume  the  buried  layer  is  approximated  by  the  same 

type  impurity  profile  used  to  approximate  the  dov/n  emitter 

region.  Therefore,  in  combination,  we  obtain  the  following 

2 

approximating  equation  for  the  intrinsic  region  of  an  I L 
type  structure : 


C(x)  = 


C ,erfc 
ol 


fe]- 


^o2®^P 


fK+xl 


[2L2j 


2^  J 


+ C _erfc 
o3 


W . -K-xl 
epi 


2L- 


+ C 
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In  Eq.  3,  , L2,  and  represent  the  impurity  atom  diffu- 
sion length  (L  = /Dt)  for  each  species;  and  is  the  epi- 

taxial layer  thickness.  Throughout  the  extrinsic  region  of 
this  structure  we  introduce  an  obvious  modification  of  Eq.  3 
we  eliminate  the  down  emitter  diffusion  profile  term: 


C(x)  = - 


Co2exp 


K+x' 

2 

W ,-K-x 

+ C ,erfc 
o3 

epi 

2L2 

L ^^3  J 

+ C. 


(4) 


Thus,  from  Eqs.  3 and  4 we  have  the  basis  of  our  phys- 
2 

ical  model  for  I L operation.  Using  these  mathematical 
relations,  an  iterative  procedure  is  used  to  establish  the 
junction  locations  (where  C (x)  = 0)  and  the  buried  layer 
edge  is  defined  by  that  location  where 


C _erfc 
o3 


W .-K-x 
epi 


2L 


= C 


3 J 


B 


(5) 


Figs.  4 and  5 illustrate  typical  impurity  atom  profiles 

2 

used  in  this  mathematical  model  for  an  I L type  semiconduc- 
tor device. 
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For  convenience.  Fig.  6 illustrates  the  notation  used 

in  this  report  to  identify  the  junctions,  and  their  respec- 

2 

tive  space-charge  layer  edges,  within  a typical  I L vertical 
transistor. 


2 . 2 Space-Charge  Layer  Characteristics — Equilibrium 

The  depletion  approximation  represents  a traditionally 
important  aspect  of  p-n  junction  theory.  This  approximation 
technique  was  first  used  in  the  analysis  of  abrupt  p-n  junc- 
tions [5] , and  its  application  continued  for  many  years. 

It  is  now  realized  that  the  depletion  approximation  has  only 
limited  applicability  for  the  abrupt  p-n  junction  [6-71,  al- 
though this  approximation  is  generally  applicable  to  linearly 
graded  p-n  junctions  [8],  and  to  diffused  p-n  junctions  [9]. 
Here  we  model  the  space-charge  regions  of  diffused  junctions 
based  upon  traditional  concepts  of  the  depletion  approximation. 
At  equilibrium,  we  assume  the  space-charge  regions  of 

2 

I L junctions  can  be  represented  by  solutions  of  Poisson's 
equation : 


d III  ^ 

^ 2 
dx 


C (x)  , 


[5]  W.  Shockley,  Bell  Syst.  Tech.  Jour.,  28,  435  (1949). 

[6]  D.P.  Kennedy,  IEEE  Trans.,  ED-22,  11  TT975). 

[7]  D.P.  Kennedy,  Solid  State  Electronics,  20 , 311  (1977). 

[8]  D.P.  Kennedy  and  R.R.  O'Brien,  IBM  J.  of  Res.,  11, 

252  (1967). 

[9]  D.P.  Kennedy,  Math.  Invest,  of  Semicond.  Devices, 

Final  Report  AF  19 (628) -5072 , AFCRL-66-358 (11) , 1967. 
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n-type  p-type 


^ni  ^pi 


n-type 


^P2  ^j2  ^n2  ^J3 


Xji  = Down  Emitter 
Xj2  = Up  Emitter 


Fig.  6 Illustration  of  the  Notation  Used  to  Identify  Each 

Junction  (and  its  respective  space-charge  layer  edges) 
in  the  Intrinsic  Region  of  an  I^L  Vertical  Transistor. 
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where  C(x)  for  intrinsic  regions  is  given  by  Eq.  3 and 
for  extrinsic  regions  by  Eq.  4.  Further,  from  Eq.  6 the 
space-charge  layer  electric  field  distribution  is  given  by 


K (x^ ,x) 


C { t)  dt  , 


(7) 


In  addition,  a second  integration  of  Eq.  7 yields  the 
electrostatic  potential  distribution  throughout  a p-n  junc- 
tion , 


\p  (x) 


C(t)dtdn  . 


(8) 


It  can  be  shown  [10]  that  both  Eq.  7 and  Eq.  8 are  expres- 
sible in  closed  analytical  form  and,  therefore,  are  easily 
calculated . 

In  addition,  at  equilibrium  it  is  reasonable  to  assume 
all  mobile  carrier  densities  within  the  p-n  junction  are 
adequately  approximated  by  a Boltzmann  type  of  distribution. 


_ 

a)  p = C^e 

b)  n = 


(9) 


From  Eq.  9a  and  Eq.  9b,  the  total  voltage  across  a diffused 
p-n  junction  space-charge  layer  at  equilibrium  is  given  by: 


[10]  M.  Abramowitz  and  J.A.  Stegum,  Handbook  of  Mathematical 
Functions,  U.S.  Deot.  of  Commerce , 1964 . 


(10) 


V = — Jin 
eq  q 


p (x  ) n (x  1 
^ pe  ne 


n2 

1 


where 


C(x)  + /C(x)2  + 4n? 


n(x)  = p(x)  = 


(11) 


In  Eq.  10,  Xp^  and  x^^  designate  the  space-charge  layer  edges 
in  the  n-type  and  p-type  material,  respectively,  at  equilibrium. 

Thus,  in  combination,  Eqs . 7,  8,  and  10  yield  the  equi- 
librium conditions  for  a diffused  p-n  junction: 


X 


a)  0 = 


Ke 


o X 


pe 


C(t)dtdn  - 


kT  . 
- ^n 


ne  ne 


p (x  ) n (x  ) 
^ pe  ne 


X 


n . 

1 

(12) 


b)  0 = 


_a_ 

K£ 

O 


pe 


ne 


C(t)dt  . 


This  mathematical  model  has  been  used  to  calculate  the 

equilibrium  space-charge  layer  edges  for  both  the  up  and  down 
2 

emitters  of  an  I L type  structure. 

Using  this  model,  we  locate  in  a given  device  the  space- 
charge  layer  edges  of  both  p-n  junctions,  and  their  respec- 
tive "built-in"  or  equilibrium  diffusion  voltages. 


V 

eq 


kT 

q 


Jin 


p (x  ) n (x  )' 
pe  ne 


2 
n . 
1 


(13) 


From  a knowledge  of  these  space-charge  layer  edge  locations, 
(x^^,Xp^)  we  also  calculate  the  equilibrium  space-charge 
layer  capacitance  of  each  junction. 
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1 


In  addition,  a knowledge  of  these  locations  provides  a 

means  to  readily  establish  the  equilibrium  electrical 

2 

base  regxon  width  within  an  I L transistor, 


f 

W = X -X 

e pe,  pe,, 


2 . 3 Space-Charge  Layer  Characteristics--Reverse  Bias 

Equation  7 and  Eq.  8 are  also  used  to  calculate  the 
reverse  biased  characteristics  of  a diffused  p-n  junction. 
With  a knowledge  of  the  "built-in"  equilibrium  voltage  (Sec. 
2.2),  the  reverse  bias  space-charge  layer  edges  are  located 
by  satisfying  the  relations 


a)  0 = C(t)dtdn 

n n 


V - V 
eq  a 


b)  0 = 


—2—  C (t) dt  . 

ice 


In  Eq.  16,  represents  the  equilibrium  junction  barrier 
potential  and  is  an  applied  reverse  biasing  voltage. 

As  used  for  the  equilibrium  case  'Sec.  2.2),  the  cal- 
culation establishes  the  space-charge  layer  edges 
for  each  p-n  junction.  Similarly,  after  calculating  these 
locations  we  establish  the  space-charge  layer  capacitances 
for  each  individual  junction. 
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(17) 


C (V. 


Ke 


X -X 

n p 


and  the  nonequilibrium  electrical  base  width  for  this  semi- 
conductor structure 


W = 
e 


X -X 

P2  Pi 


(18) 


2 . 4 Reach-Through  Voltage  Calculations 

Another  application  of  Eqs.  7 and  8 yields  the  reach- 

through  voltages  for  this  semiconductor  device.  Specifically, 

the  reverse  biasing  voltage  required  to  produce  space-charge 

layer  reach- through  within  the  transistor  base  region.  It 

has  been  shown  [2]  that  this  particular  electrical  parameter 

offers  many  difficulties  with  respect  to  bipolar  transistor 

fabrication  reproducibility,  particularly  in  narrow-base 

configurations.  It  will  be  shown  here  that  space-charge 

layer  reach-through  can  be  a similar  problem  in  the  fabrica- 
2 

tion  of  I L structures. 

During  up  operation,  the  reverse  biased  up  collector 
space-charge  layer  edge  (X  in  Fig.  6)  can  reach-through 

Pi 

to  the  up  emitter  space-charge  layer  edge  (X  in  Fig.  6) . 

P2 

In  this  type  of  situation,  the  total  up  collector  space- 
charge  layer  can  be  established  using  Eq.  7, 


0 = 


X 

^ P2 


n. 


C(x)dX  . 


(19) 
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Thereafter,  using  Eq.  8,  we  establish  the  up  collector 
punch- through  voltage. 


RT  Kc 


P2  n 


C (t) dtdn 


o ^'X 

ni  n^ 


Substantial  changes  have  been  found  in  the  location 

of  X for  forward  bias  and  equilibrium  operation  of  the 

up  emitter.  As  a consequence,  the  up  collector  reach-through 

voltage,  V , changes  for  these  two  modes  of  device  opera- 
RT 

tion,  and  it  is  necessary  to  specify  two  different  reach- 

through  voltages:  emitter  is  at  zero  bias 

and  when  the  up  emitter  is  operational. 

Assuming  the  up  emitter  junction  space-charge  layer 

edges  are  at  their  equilibrium  locations  (X  and  X ) , 

^ pe2  ne2 

we  satisfy  the  relation: 


0 = C(t)dt 

X 


Thereafter,  from  Eq.  8 we  have  the  reach-through  voltage. 


^2  .n 

V = -3- 

"l  "l 


c(t)dtdn  . 


When  the  up  emitter  junction  is  forward  biased  its 
space-charge  layer  width  is  exceedingly  narrow.  For  pur- 
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poses  of  this  punch-through  calculation,  little  error 


arises  by  assuming  X = X.  and  therefore, 

^2  ^2 


X . 

^2 

0 = C(t)dt  , 

''X 


C ( t) dtdn 


X X 
"l  "l 


Contrasting  with  up  emitter  operation,  little  change 


arises  in  the  down  emitter  space-charge  between  equilibrium 


and  forward  bias.  This  situation  arises  from  the  large 
impurity  atom  gradient  associated  with  this  p-n  junction. 
For  this  reason,  it  becomes  necessary  to  define  only  one 
reach-through  voltage  for  down  transistor  operation, 
and  this  voltage  is  approximated  by  the  equations 


C (x)  dx 


in  combination  with 


1 rh 


f -L 
X 


C (t) dtdn 


X J X 
*^2  '^2 
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2 . 5 Junction  Avalanche  Breakdown 

The  avalanche  breakdown  of  a diffused  p-n  junction  can 
be  readily  modeled  for  any  arbitrary  space-charge  layer  im- 
purity atom  distribution  [11-13].  This  calculation  is 
accomplished  by  satisfying  the  relation 


^n 


where,  from  Eq.  8, 


C(n)  = ^ C(t)dt 
o X 

n 


0 = C(t)dt  . 

■'X 


After  establishing  the  space-charge  layer  edges  satisfying 
Eqs.  27-29,  the  avalanche  breakdown  voltage  is  given  by 


P rn 


C (x) dxdn 


'o  'X  'X 
n n 


It  will  be  shown  that  the  "down"  operation  collector  junction 
will  scmetirrEs  undergo  avalanche  breakdown  before  reach-through 
In  this  situation,  measurement  of  the  collector 

[11]  D.P.  Kennedy  and  R.R.  O'Brien,  IRE  Trans.,  ED-9,  478 
(1962) . 

[12]  D.P.  Kennedy  and  R.R.  O'Brian,  IBM  Jour. , 422  (1965). 

[13]  D.P.  Kennedy  and  R.R.  O'Brien,  IBM  Jour. , 10,  2]3  (1966). 
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breakdown  voltage  will  be  shown  to  offer  an  experimental 

means  for  gaining  valuable  information  concerning  physical 

2 

properties  of  an  I L device. 


2.6  Ooen-Base  Breakdown  BV 


As  yet,  open-base  breakdown  calculations  have  not 

2 

been  implemented  in  our  mathematical  model  for  I L opera- 
tion. Nevertheless,  this  has  been  previously  accomplished 

for  conventional  bipolar  transistor  operation  [14],  and 

2 

these  techniques  are  directly  applicable  to  the  I L type 
structure . 

Specifically,  open  base  breakdown  occurs  when  aM  = 1 
(hence  B = °°)  , where  a is  the  grounded  base  current  gain, 
and  M is  the  current  multiplication  occurring  within  a 
collector  junction  space-charge  layer.  Thus,  open-base 
breakdown  can  be  calculated  by  satisfying  the  relation. 


(1-a)  = 


dn  , 


where  C(n)  is  given  by  Eq.  28  and,  in  addition,  the  relation 


0 = C(t)dt  . 

X 


[14]  D.P.  Kennedy  and  R.R.  O'Brien,  Int.  Jour,  of  Elect., 
18,  133  (1965). 


Eqs.  31  and  32  establish  the  junction  space-charge  layer 
edges  at  open-base  breakdown.  After  completing  this  cal- 
culation, the  breakdown  voltage  is  given  by 


BV 


CEO 


X 

r P 


C (t) dtdn  . 


(33) 


As  previously  stated,  this  model  has  not,  as  yet,  been 

2 

implemented  in  our  analysis  of  the  I L type  structure.  This 

is  not  intended  to  imply  that  we  feel  such  a calculation  is 

unimportant.  During  "up"  operation  substantial  collector 

multiplication  could  occur;  even  at  the  low  operating  voltages 

of  this  semiconductor  device.  This  calculation  is  suggested 

because  a recent  publication  [15]  indicated  that  open-base 
2 

breakdown  in  I L devices  is  a problem  of  concern  for  device 
designers. 


2 . 7 and  Stored  Intrinsic  Base  Charge 

It  has  been  shown  that  the  minority  carrier  dis- 

tribution within  an  intrinsic  base  layer  can  be  approximated 
by  the  relation: 


n (x) 


J 


n 

qDn^(x) 


C (x) dx  . 


(34) 


For  greater  accuracy,  we  have  introduced  into  Eq.  34  the 
average  base  region  electron  diffusion  coefficient,  where 

[15]  J.L.  Saltich,  W.E.  George,  and  J.G.  Soderberg,  IEEE 
Trans. , SC-11,  478  (1976). 
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kT 


(35) 


D 

n 


q 


M (h)dh  . 
X ^ 


From  these  expressions,  the  required  to  attain  a given 

emitter  junction  current  density  is  given  by 


V 


kT 


EB 


in-i 


n 


— 2 
qD  n 
^ n 1 


C(x)dx  + 1 


(36) 


X 


Our  implementation  of  this  model  was  based  upon  an 

assumption  that  I L "up"  operation  does  not  produce  high 

minority  carrier  injection  levels  into  the  transistor  base 

region.  After  studying  the  results  of  several  calculations 

it  was  found  that  this  assumption  is  sometimes  unjustified; 

further  details  of  this  study  are  presented  in  Section  3. 

of  the  present  report.  As  a consequence,  a new  model  is 

being  developed  to  replace  Eq.  34  for  V__  calculations. 

Despite  this  deficiency  in  our  model,  it  is  shown  (Section 

2 

3.3)  that  other  concerns  of  I L analysis  and  design  are  far 

more  important  at  the  present  time. 

In  addition  to  calculating  V„_  (Eq.  36)  this  model  also 

BB 

yields  the  stored  charge  due  to  minority  carriers  within  the 
intrinsic  base  region: 


r ^2 


^int  D 

n 


rn 


C(n) 


C(t)dtdn  . 


(37) 
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Further,  from  p.  .we  also  obtain  the  intrinsic  base 
int 

region  recombination  current  J . . 

ri 


2 . 8 Extrinsic  Base  and  Up  Emitter  Region  Stored  Charge 

2 

We  next  model  a property  of  the  I L transistor  that 
is  important  to  its  switching  speed:  the  stored  extrinsic 

charge  during  "up"  operation.  Here  we  use  relatively  con- 
ventional assumptions  concerning  the  associated  mechanisms. 
It  is  assumed  the  extrinsic  base  region  minority  carrier 
current  is  attributable  to  the  combined  mechanisms  of 
drift  and  diffusion, 


J = qD  ^ + qp  nf (x)  , 
n ^ n dx  ^^n  ' 


and  that  continuity  of  this  injected  electron  current  is 
governed  by 


where  represents  the  base  region  minority  carrier  life- 
time. This  condition  of  continuity  is  recognized  as  the  low 
injection  limit  of  conventional  SRH  recombination  processes. 

After  combining  Eqs.  38  and  39  we  obtain  the  differen- 
tial equation 


dn.  1 dC  n(x)_,^ 
dx2  C(x)  dx  - -^2  = 0 ' 

n 
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I 
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1 

t 


1 


1 

[■ 


where  the  electric  field  ^(x)  in  Eq.  38  is  obtained  from 
the  quasi-neutral  approximation 


C(x) 


M 1 dC (x) 
q C(x)  dx 


(41) 


In  the  past,  two  different  approaches  have  been  used 
to  solve  Eq.  40:  first,  an  infinite  series  type  solution 

[16]  and,  second,  a finite  difference  type  of  solution. 

The  first  method  requires  a little  more  computer  time,  yet 
it  does  not  exhibit  the  errors  known  to  be  encountered  in 
any  finite-difference  type  solution.  The  second  method  is 
computationally  fast,  it  is  easier  to  implement,  although 
this  technique  can  lead  to  round-off  error,  and  instabilities, 
that  sometimes  produce  problems.  In  our  implementation  of 
this  model  we  adopted  the  finite-difference  approach.  It 
is  planned  at  a later  time  to  implement  an  infinite  series 
method  for  solving  Eq.  40  and,  thereafter,  draw  comparisons 
between  these  two  computational  techniques. 

Conventional  methods  have  been  employed  to  obtain  this 
solution.  Equation  40  is  first  written  in  its  finite- 
difference  form,  yielding  the  well-known  tri-diagonal  matrix. 
Thereafter,  this  system  of  equations  is  solved  using  Gaussian 
elimination  [17],  and  backward  calculation  of  the  resulting 
minority  carrier  (electron)  distribution. 


[16] 

[17] 


D.P.  Kennedy  and  P.C.  Mur ley,  IRE  Trans 
C.E.  Pearson,  Handbook  of  Applied  Math, 
Reinhold,  New  York  (1974 ) . 


ED  9,  136  (1962)  . 
Van  Nostrand 
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I 
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Boundary  conditions  used  for  our  solution  of  Eq.  40 
are  based  upon  an  applied  at  the  up  emitter  junction, 

and  an  infinite  surface  recombination  velocity  at  the  semi- 
conductor-oxide interface.  From  previous  calculations  of 
this  type  [17],  during  "up"  operation  it  was  suspected  that 
the  "built-in"  base  region  electric  field  was  sufficient 
to  keep  an  appreciable  number  of  minority  carriers  from 
this  interface  region.  For  this  reason,  by  assuming  an 
infinite  surface  recombination  velocity  we  are  taking  a 
most  pessimistic  view  concerning  the  influence  of  the  surface 
upon  device  operation.  j 

An  identical  approach  is  used  to  calculate  stored  ] 

] 

charge  within  the  "up"  emitter  region  of  this  structure.  i 

■I 

In  our  model  we  assume  an  infinite  recombination  velocity  | 

at  the  substrate-epitaxy  interface.  Calculations  indicate  I 

that  the  quasi-neutral  electric  field  produced  by  the  buried  | 

layer  negates  any  assumption  concerning  this  interface  region.  " | 

Calculations  indicate  that  near  this  interface  we  obtain  j 

minority  carrier  density  that  is  10^^  to  10^^  below  that  ] 

I 

found  near  the  "up"  emitter  junction.  | 

After  solving  Eq.  40  and  obtaining  a predicted  minority  i 

j 

carrier  distribution,  the  base  region  stored  charge  is  ob-  J 

\ 

tained  from  i 


^ext 


q 


n (x) dx 


(42) 


for  the  extrinsic  base  layer,  and  from  | 
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Pe  = q 


w 

r epi 

X 


p (x)  dx 


(43] 


for  the  emitter  region  under  "up"  operation. 


2 

3 . 0 Analysis  of  the  I L Vertical  Transistor 

Using  the  foregoing  model,  many  important  electrical 
properties  of  an  I L type  transistor  can  be  readily  es- 
tablished. Resulting  from  this  study  is  a better  under- 
2 

standing  of  t l operation,  some  insight  into  basic  dif- 
ferences between  this  type  device  and  conventional  bipolar 
transistor  operation,  and  an  understanding  of  many  asso- 
ciated problems  relating  to  device  fabrication.  Thereby, 
a view  is  obtained  of  the  capabilities  inherent  in  this 
proposed  modeling  technique,  and  its  potential  value  for 
modeling  the  influence  of  radiation  induced  mechanisms  upon 
the  electrical  properties  of  an  I L type  semiconductor 
device . 

2 

Figure  7 illustrates  the  I L structure  adopted  for 

the  present  study;  this  structure  is  similar  to  that  used  in 

2 

a previous  investigation  of  I L operation  [18].  The  changes 
introduced  here  resulted  from  an  initial  evaluation  of  the 
previously  used  model;  calculations  indicated  emitter  to 
collector  reach-through  at  equilibrium,  a situation  considered 

2 

[18]  J.A.  Niehaus,  Computer  Aided  Analysis  of  I L,  Univ.  of 
111.  Report  UlLU-ENG  75-2224,  Joint  Services  Contract 
DAAB-07-72-C-0259,  Sept.  1975. 
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impractical.  For  purposes  of  illustration  we  have  adopted 

3 

the  following  impurity  atom  concentrations  (atom/cm  ) : 


TABLE  I 

Fabrication  Process  Parameters  Assumed  for  the  Design-Center 


Device 

23 

Si  = 

19 

Cq2  - 5x10  ^ 
S3  = 

S = 10^^ 


It  is  again  emphasized  that  the  exceedingly  large  concentra- 
tion used  here  for  is  not  intended  to  represent  a 

realistic  value.  Instead,  this  large  concentration  is  used 
to  approximate  the  impurity  atom  distribution  obtained  in 
the  active  regions  of  devices,  due  to  nonlinear  diffusion 
mechanisms . 


3. 1 Equilibrium  Calculations 

At  equilibrium  we  calculate  the  following  physical  and 

2 

electrical  characteristics  for  this  illustrative  I L tran- 


sistor; 
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TABLE  II 


Equilibrium  Characteristics  for  the  Design-Center  Device 


Junction  #1 
X = 1.14 9 urn 

nei 

X.  = l,200ym 
3l 

X = 1.297ym 


= 0.834  volts 


C^  = 0.014  Pfd 


Junction  #2 


X = 1.440ym 

pe2 

X ■ =1. SOOym 


X = 1. 960um 

ne2 

V = 0.695  volts 

eq2 

C^  = 0.239  Pfd 


Transistor 


Wp  = 0.400ym 

W = 0 . 14  3ym 
e 


An  observation  from  this  calculation  is  the  difference 
between  the  physical  base  width  (W^)  and  the  electrical 
base  width  (W^) . In  this  particular  structure  only  about 
35%  of  the  physical  base  region  is  charge  neutral.  This 
situation  could  produce  a substantial  fabrication  repro- 
ducibility problem.  For  this  reason,  additional  equilibrium 
calculations  provide  insight  into  the  factors  influencing 


these  particular  device  parameters. 
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X = 1.423ym 

pe  ^ 

X.  = 1.600ym 
^2 

X = 1.950um 

ne2 

V = 0.695  volts 

^2 

Co  = 0.236  Pfd 


Transistor 


W = 0.400ym 
P 

Wg  = 0.098ym 


This  proposed  modification  produces  only  a small 
reduction  of  the  junction  space-charge  layer  capacitances 
(due  to  an  increase  of  space-charge  layer  width) , but  it 
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does  produce  a large  change  of  electrical  base  width. 

19  3 

Specifically,  reducing  from  5x10  (atoms/cm  ) to 

19  3 

10  (atoms/cm  ) reduces  the  equilibrium  electrical  base 
region  from  .143um  to  0.098ijm;  a reduction  of  about  31%. 

It  has  been  shown  [19]  that  in  an  "up"  direction  of 
operation,  emitter  region  back  injection  represents  an 
important  source  of  difficulty.  In  addition,  it  was  shown 
that  a buried  layer  can  be  used  to  minimize  operating 
difficulties  in  this  region;  this  layer  reduces  minority 
carrier  recombination  at  the  substrate-epitaxy  interface. 

As  a consequence,  there  is  a tendency  for  I L engineers 
to  move  this  buried  layer  very  close  to  the  adjacent  p-n 
junction;  this  design  method  is  intended  to  improve  emitter 
injection  efficiency  in  the  "up"  direction.  For  this  reason, 
the  consequences  of  such  a design  modification  are  esta- 
blished for  this  device  at  equilibrium. 

Calculations  were  performed  assuming  the  previously 
stated  design-center  fabrication  parameters  of  this  illus- 
trative device.  Table  I,  except  the  buried  layer  location 
was  changed  to  l.SOym.  This  design  modification  was  accom- 
plished by  assuming  in  our  model  a reduction  of  epitaxial 
layer  thickness — from  3.0um  to  2.2um.  Such  a change  would 
locate  the  -turied  layer  close  to  the  emitter  junction  space- 
charge  layer  edge  during  "up"  operation. 


[19]  F.M.  Klaassen,  IEEE  Trans. , ED-22,  145  (1975). 


Calculations  show  that  this  change  produced  an  increase 
of  equilibrium  space-charge  layer  capacitance  (from  0.244  Pfd 
to  0.36  Pfd),  with  little  significant  change  of  electrical 
base  width.  Hence,  from  a device  design  point  of  view,  a 
reduction  of  epitaxial  layer  thickness  (from  3.0um  to  2.2ym) 
appears  to  be  a reasonable  design  modification. 

In  a previous  study  [2]  it  was  shown  that  epitaxial 
layer  thickness  represents  an  important  source  of  fabrica- 
tion process  variability.  From  wafer  to  wafer  an  observable 
difference  is  found  in  the  epitaxial  thickness  produced  by 
a given  reactor.  Similarly,  changes  of  epitaxial  layer 
thickness  can  also  be  observed  across  any  given  slice  of 
silicon  produced  in  a given  reactor.  In  this  previous  study 
a ±10%  thickness  change  could  be  observed  in  measurements 
taken  upon  a large  number  of  epitaxial  layers.  For  this 

reason,  it  is  important  to  consider  the  consequences  of 

2 

inadvertent  epitaxial  thickness  changes  upon  this  I L device, 
when  it  is  assumed  to  be  fabricated  on  a 2.2ym  epitaxial 
layer . 

Calculations  were  performed  assuming  the  buried  layer 

starts  at  X.  (1.6y)  and  the  epitaxial  layer  is  2.0um  thick. 
^2 

Following  are  the  results  of  these  calculations. 


Junction  #1 


X 

ne^ 

= 1.449ym 

X . 

= 1.200ym 

^1 

X 

= 1.297ym 

pei 

V 

= 0.834ym 

eqi 

^1 

= 0.014  Pfd 

X 

pe 


V 

eq 


Junction  #2 
= 1.466ym 

2 

= l,600ym 
= 1.621ym 

2 

= 0 . 7 32yiTi 

2 


C2  = 0.799  Pfd 


Transistor 
Wp  = 0.400ym 
= 0,169ym 

Table  IV  shows  the  reduction  of  epitaxial  layer 
thickness  producios  an  increase  of  impurity  atom  gradient 
and,  as  a consequence,  a reduction  of  the  equilibrium 
space-charge  layer  width  (X  -X  ).  Therefore,  this 
modification  of  the  epitaxial  layer  thickness  produces  an 


increase  of  electrical  base  width  (from  0.143um  to  0.169ym) 


Junction  #1 


X = 1.135ijm 

X . = 1 . 200pm 

^1 

X = 1.385pm 

Pi 

V = 1.5  volts 
CB 

= 0.008  Pfd 


Junction  #2 
X = 1.440pm 

P2 

X.  = 1.600pm 
^2 

X = 1 . 960pm 
VgE  = 0.000  volts 
C2  = 0.239  Pfd 


Transistor 

W = 0.400pm 
P 

Wg  = .056pm 


It  should  be  observed  that  this  calculation  is  based 
upon  an  "up"  emitter  that  is  at  equilibrium  = 0) . 

The  biasing  voltage  assumed  on  this  "up"  collector 
junction  (V„q  = 1.5V)  produces  an  electrical  base  width 


that  may  be  too  small.  One  consequence  of  this  situation 
is  the  overall  reproducibility  problem  in  the  diffusion 
process  [20] . For  this  reason,  similar  calculations  have 
been  made  for  this  device  assuming  an  "up"  collector  bias 
of  1.0  volts.  Following  are  the  results  of  such  calcula- 


tions : 


TABLE  VI 


Calculated  Characteristics  for  the  Design-Center  Device 
Assuming  = l.OV;  = 0 

Junction  #1 

X = 1.138ym 
"l 

X . = 1 . 200ym 

=>1 

X = 1.355ym 


= 1.000  volts 
C^  = 0.009  Pfd 


Junction  #2 


X = 1.440ym 
^2 

X.  = l.eOOym 

32 

X = 1.960ym 


^EB  = ^ 

= 0.239  Pfd 


Transistor 


V'J  = 0.400ym 
P 

= 0.086ym 


[20]  D.P.  Kennedy,  IBM  Journ. , 5,  331  (1961). 
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From  these  calculations,  a reduction  of  0.5  volts 

bias  on  the  up  collector  produces  an  increase  of  electrical 

base  width  (assuming  = 0)  from  0.056ym  to  0.086um. 

bo 

In  a practical  device  fabrication  situation  it  is  indeed 
questionable  that  such  a device  would  represent  a prac- 
tical design.  Small  inadvertent  variations  of  boron  im- 
purity atom  surface  concentration  could  produce  a condi- 
tion of  reach-through  in  this  mode  of  operation. 


3.3  Reverse  Biased  Up  Collector--V„„?^0 
■ " ' — ^Eo 

Assuming  an  up  collector  junction  bias  of  1.0  volts, 

we  next  calculate  the  required  to  attain  a specified 

bo 

collector  current. 


TABLE  VII 


Calculations  for  the  Design-Center  Device 

— K K - ■ - ■ — 


I^/collector 

J 

-c 

V 

-EB 

( M-amps) 

(amps/cm  ) 

volts 

5.0 

25.0 

0.603 

10.0 

50.0 

0.621 

50.0 

250.0 

0.663 

100.0 

500.0 

0.681 

200.0 

1000.0 

0.699 

Table  VII  illustrates  a situation 

encountered 

in 

I^L 

transistors  that 

is  not  characteristic 

of 

traditional 

bi- 

polar  operation. 

The  small  impurity  atom 

gradient 

at 

the 

up  direction  emitter  places  a limitation 

upon  the  i 

minority 
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current  density  obtainable  from  this  junction.  This  limita- 
tion arises  from  the  small  equilibrium  voltage  (V  = 0.695 
volts  in  our  illustrative  example) . In  recent  publications 
[18-19]  it  was  assumed  that  the  up  direction  emitter  junc- 
tion  would  yield  a current  density  of  1000  Amps/cm  with  an 
applied  of  about  0.750  volts;  these  assumptions  would 

be  t' oical  of  the  down  direction  emitter  junction.  In  the 
up  direction  a similar  current  density  might  not  be  at- 
tained without  specific  design  modifications. 

In  addition  to  the  equilibrium  junction  voltage,  another 
factor  must  be  taken  into  consideration:  up  emitter  junc- 

tion space-charge  layer  contraction,  due  to  an  applied 
forward  bias.  A substantial  increase  of  electrical  base 

width  results  from  the  V applied  to  the  up  emitter  junction 

EB 

and,  as  a consequence,  the  up  electrical  base  width  is  substan- 
tially wider  than  the  down  electrical  base  width. 

Calculations  show  that  by  moving  the  buried  layer 
slightly  into  the  equilibrium  space-charge  layer  of  the  up 
emitter  junction,  an  increase  is  obtained  in  the  built-in 
voltage.  For  example,  from  TABLE  II,  the  equilibrium  up 
emitter  junction  terminates  at  X = 1.96)jm.  A design 

"2 

modification  placing  X.  = 1.7um  (slightly  inside  the  equi- 

ps 

librium  space-charge  layer) , the  equilibrium  junction  vol- 
tage is  increased  from  0.695  volts  to  0.768  volts. 
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Electrical  Base  Width  (V„  =1.0  volt;  = 100  p-amp/ 
Collector) 


It  was  previously  stated  that  a substantial  change  of 
electrical  base  width  results  from  forward  biasing  the 


emitter  junction.  Assuming  a collector  current  of  100  p-amps, 


the  following  physical  and  electrical  characteristics  were 
calculated  for  this  proposed  design-center  device  (TABLE  I) : 


TABLE  VIII 


Calculated  Characteristics  Assuming  = 1.0  V;  = lOOpa 


Junction  #1 


X 

'^l 

1.138pm 

X . 

n 

1. 200pm 

X 

Pi 

1 . 355pm 

II 

m 

u 

> 

1.00  volts 

^1  = 

0.009  Pfd 

Junction  #2 

X 

P2 

1 . 596pm 

X . = 

=>2 

1 . 600pm 

X 

■^2 

1 . 610pm 

'^EB 

0.681  volts 

^2  = 

85.7  Pfd 

Transistor 

W 

P 

0 . 400pm 

W 

e 

0 . 241pm 

f 


I 
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From  TABLE  VIII,  when  this  I L device  is  operating  in 

the  up  direction  (J^  = lOOpa)  we  have  an  electrical  base 

width  of  0.298Mm,  which  is  larger  than  encountered  in 

traditional  bipolar  operation. 

It  should  also  be  noted  that  this  up  emitter  junction 

has  a very  large  space-charge  layer  capacitance  at  this 

particular  biasing  voltage.  It  is  suggested  that  this 

capacitance  could  have  a significant  influence  upon  the 

2 

switching  speed  of  an  I L structure. 


i 


3. 5 Intrinsic  Base  Region  Stored  Charge 

Using  this  model,  the  intrinsic  base  region  stored 
charge  has  been  calculated  throughout  a wide  range  of 
collector  currents.  As  expected,  this  stored  charge  is 
linearly  dependent  upon  the  collector  current: 

TABLE  IX 

Intrinsic  Base  Region  Stored  Charge  (Vgg  =1.0  volts) 


I^/Collector 

^c 

Pj^^^/Col  lector 

(y-amps) 

2 

(amps /cm  ) 

(Pico-Coulombs ) 

5.00 

25.0 

3.1x10"^ 

10.0 

50.0 

6.2x10“^ 

50.0 

250 

3.1x10"'^ 

100 

500 

6.2x10”'^ 

200 

1000 

1.2x10“^ 

Included  in  this  model  is  the  capability  of  calculating 
the  minority  carrier  distribution  between  the  emitter  and 
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collector  junctions.  Figure  8 illustrates  this  calculated 
minority  carrier  distribution  assuming  a collector  current 
of  100  M-amps. 

Figure  8 provides  substantial  insight  into  the  in- 
fluence of  the  "built-in"  base  region  electric  field.  It 
has  been  suggested  [21]  that  during  up  operation  the 
"built-in"  electric  field  in  this  device  is  directed  to 
oppose  carrier  transport  from  the  emitter  to  the  collector. 
As  a consequence,  this  "built-in"  field  is  believed  to 

produce  excessive  base  region  charge  storage  and,  there- 

2 

fore,  reduce  the  switching  speed  of  I L structures. 

Figure  8 indicates  that  the  stored  charge  in  this  intrinsic 
base  region  differs  little  from  that  of  a homogeneous  base 
device.  As  a consequence,  it  is  difficult  to  understand 
how  this  particular  stored  charge  mechanism  could  con- 
tribute to  the  switching  speed  problems  associated  with  the 
2 

I L technology. 

3 . 6 Extrinsic  Base  Region  Stored  Charge 

In  the  device  we  have  selected  for  analysis,  the 
"built-in"  extrinsic  base  region  electric  field  appears  to 
be  of  significant  advantage.  In  particular,  this  electric 
field  tends  to  reduce  the  extrinsic  base  region  stored 
charge  and,  thereby,  reduces  the  base  region  recombination 
current.  As  a consequence,  the  extrinsic  base  region 

[21]  C.  Thornton,  ECOM  (private  communication). 
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minority  carriers  accumulate  near  the  up  emitter  junction 
face,  and  the  density  of  these  minority  carriers  decreases 
rapidly  with  an  increase  of  distance  from  this  location. 

Fig.  9. 

Calculations  similar  to  Fig.  9 have  been  undertaken 
throughout  a wide  range  of  minority  carrier  lifetime 
(10~®<t<10  with  little  change  in  this  distribution. 

From  these  calculations,  and  the  distances  involved  rela- 
tive to  the  minority  carrier  diffusion  length,  it  is 
evident  that  the  distribution  in  Fig.  9 is  not  necessarily 
attributable  to  recombination.  This  distribution  is  assumed 
a consequence  of  the  "built-in"  electric  field  produced  by 
the  base-region  impurity  atom  distribution. 

From  these  calculations  we  have  also  established  the 
extrinsic  base  region  stored  charge  throughout  the  operating 
range  of  this  semiconductor  device: 


TABLE  X 


Extrinsic  Base  Region 

Stored  Charge 

= 1.0  Volts) 

J /Collector 

J 

— D 

P 

— 

— c 

^xt 

(M-amps) 

o 

(amps/cm  ) 

(Pico-Coulombs ) 

5.0 

25.0 

2.07x10“^ 

10.0 

50.0 

4.15x10'^ 

50.0 

250.0 

2.07x10'^ 

100.0 

500.0 

4.15x10"^ 

200.0 

1000.0 

8. 314x10”^ 

i 

f!; 
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Comparing  this  extrinsic  stored  charge  with  the  total 
intrinsic  stored  charge  (TABLE  IX) , it  is  found  that  a 
substantially  greater  part  of  the  total  base  region  stored 
charge  is  in  the  extrinsic  region  of  this  structure.  A 
similarly  greater  part  of  the  base  region  recombination 
current  will  also  be  attributable  to  this  extrinsic  base 
region. 


. Ji. 


3 . 7 Up  Emitter  Region  Stored  Charge 

Many  calculations  have  been  undertaken  to  establish  the 
up  emitter  region  minority  carrier  distribution,  and  the 
stored  charge  within  this  part  of  the  structure.  These 
calculations  clearly  show  a need  for  locating  the  buried 
layer  very  close  to  the  junction  space-charge  layer;  this 
conclusion  is  in  agreement  with  a study  of  this  subject  by 
Klaasen  (19].  For  purposes  of  illustration,  Fig.  10  shows 
this  calculated  emitter  region  carrier  distribution  for 
four  different  buried  layer  locations. 

It  should  also  be  stated  that  these  calculations  show 
two  speciij.c  advantages  to  be  gained  by  locating  the  buried 
layer  close  to  the  up  emitter  space-charge  layer  edge: 

1.  An  increase  of  carrier  gradient  and,  hence, 
a decrease  in  their  integrated  density. 

2.  Some  decrease  in  the  density  of  carriers 
at  the  space-charge  layer  edge,  due  to  an 
increase  of  doping  at  that  location. 
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Minority  Carrier  Distribution  in  the  Up  Emitter  Region 
(V,,„  = l.OV;  U = lOOua)  . 


In  combination,  following  are  emitter  region  stored- 


charge  calculations  for 


the  X . locations  shown 


this  semiconductor  device, 
in  Fig.  10. 


assuming 


TABLE  XI 


Up  Emitter  Region  Stored  Charge 
X . (pm) 

2.1 

2.0 

1.9 

1.8 

1.7 


(I^  = lOPMa/collector) 

p (Pico-Coulombs) 
tie— i i. 

2.57xl0"^ 

1.50x10“^ 

4.7x10"^ 

8.7x10”'* 

2.2x10”^ 


From  these  calculations  (TABLE  XI)  it  is  evident 

that  the  buried  layer  location  is  a critical  parameter 
2 

in  the  I L transistor.  All  calculations  shown  in  TABLE 
11.0  were  accomplished  by  introducing  into  this  model  a 
variation  of  epitaxial  thickness;  thereby,  we  are  im- 
plying a necessity  to  maintain  a specified  epitaxial 

2 

thickness  in  the  I L fabrication  process. 


3 . 8 Down  Operation  Reach-Through  Voltage 

In  Section  3.5  it  was  shown  that  the  buried  layer 
in  an  I L transistor  must  be  maintained  at  a location 
that  is  relatively  close  to  the  up  emitter  junction. 

This  particular  requirement  will  signif icantly  influence 
the  down  operation  reach  through  voltage.  For  this 
reason,  TABLE  XII  shows  this  calculated  voltage  through- 
out a range  of  buried  layer  locations: 
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TABLE  XII 

"Down"  Operation  Reach-Through  Voltage 


(ym) 

*3 

(volts) 

2.1 

16.3 

2.0 

14.1 

1.9 

12.6 

00 

10.9 

1.7 

9.51 

From  TABLE  XII,  there  appears  to  be  little  dif- 
ficulty associated  with  the  down  collector  reach-through 
voltage.  Throughout  the  range  of  interest  for  locating 
the  buried  layer,  this  reach-through  voltage  is  shown  to 

always  remain  above  the  normal  operating  voltage  for  this 
2 

I L structure. 

From  this  reac  i- through  calculation,  one  might  be 
tempted  to  use  measurements  of  this  parameter  as  an  indi- 
cator of  the  distance  X-  -X.  . It  can  be  shown  that  this 
voltage  is  also  dependent  upon  the  base  region  impurities 
and,  therefore,  reach-through  measurements  are  of  little 
value  for  this  particular  evaluation  technique. 


3 . 9 Up  Emitter  Junction  Avalanche  Breakdown 

In  the  particular  structure  chosen  for  analysis, 
reverse  bias  avalanche  breakdown  will  not  be  observed; 
reach  through  to  the  up  collector  occurs  at  a lower  voltage 
than  avalanche  breakdown.  Nonetheless,  avalanche  breakdown 
measurements,  in  conjunction  with  their  calculated  values, 


y 
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can  yield  important  information.  Such  measurements  are 

2 

not  undertaken  on  completed  I L devices  but,  instead,  upon 
test  devices  containing  no  up  collector  junction.  Test 
devices  of  this  type  are  usually  constructed  by  subjecting 
them  to  all  heat  cycles  used  for  completed  device  fabrica- 
tion, yet  without  introducing  the  n-type  up  collector  im- 
purities. Table  Xlll  shows  the  calculated  reverse  break- 

2 

down  voltage  for  this  type  of  I L test  structure. 


TABLE  XIII 


Avalanche  Breakdown  Voltage  for  I L Test  Structures 


X.  (pm) 
2.6 
2.1 
2.0 
1.9 
1.8 
1.7 


V (volts) 
-av 


30.2 
26.0 

23.1 

20.1 

17.3 
17.2 


These  calculations  show  that  the  avalanche  break- 
down voltage  is  exceedingly  dependent  upon  the  buried 

layer  location,  relative  to  X . . It  can  also  be  observed 

^2 

from  these  calculations  that  when  X.  • l.Svjm,  approx- 

^3 

imately,  the  structure  undergoes  breakdown  in  a manner 

similar  to  an  abrupt  junction--V^^  is  nearly  independent 

of  X . . 

^3 

Similar  breakdown  calculations  throughout  a range  of 
base  impurity  atom  density  (and  background  doping)  will 
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show  the  particular  value  of  this  type  I L test  structure. 

Specifically,  calculations  show  that  the  up  emitter 

avalanche  breakdown  voltage  is  relatively  insensitive  to 

all  fabrication  parameters  except  the  distance  X.  -x . . 

As  a consequence,  this  technique  can  be  used  to  estimate 

2 

the  distance  X . -X . in  experimental  I L devices. 

^3  ^2 

3.10  Up  Operation  Current  Gain 

Current  gain  calculations  seldom  offer  the  accuracy 

obtained  from  other  device  calculations.  Little  (or  no) 

information  is  available  concerning  the  m.inority  carrier 

lifetime  in  any  active  region  of  a bipolar  transistor. 

2 

Therefore,  for  the  present  I L analysis  a guess  must  be 
made  concerning  the  minority  carrier  lifetime,  although 
this  guess  can  be  drawn  from  experience  with  similar  cal- 
culations for  bipolar  transistor  operation.  During  bipolar 
transistor  calculations,  agreement  between  experiment  and 
theory  is  seldom  obtained  without  assuming  a minority 
carrier  lifetime  of  about  10  seconds.  For  this  reason, 

a similar  guess  has  been  made  for  all  active  regions  of 
2 

the  I L structure  under  consideration. 

In  our  calculation  of  current  gain  we  assume  all  collec- 
tor junctions  are  connected  together.  Thus,  the  total 
current  gain  for  this  device  is  given  by 
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where  is  the  electric  current/collector,  and  Ig  is  the 
total  base  current.  In  our  calculation  of  I_  we  take  into 

D 

account  the  total  base  region  recombination  current  and  the 
back  injection  current  of  the  "up"  emitter  junction.  Table 
XIV  lists  the  results  of  this  current  gain  calculation. 


TABLE  XIV 

Current  Gain  for  this  Semiconductor  Device  (I^^  = 10 Oya/col lector) 
X . (p-m)  B 


11.0 


15.6 


34.5 


67.7 


69.0 


Table  XIV  shows  the  reason  particular  emphasis  has 
been  placed  upon  the  location  of  an  I^L  buried  layer,  rela- 
tive to  the  up  emitter  junction.  These  calculations  clearly 
show  that  the  current  gain  during  up  operation  is  critically 
dependent  upon  this  buried  layer  location,  and  that  this 
location  must  be  controlled  during  device  fabrication. 

The  foregoing  calculations  are  not  intended  to  imply 

that  the  present  model  considers  all  mechanisms  influencing 

2 

current  gain  in  an  I L structure.  It  is  indeed  obvious 
that  peripheral  regions  of  the  up  emitter  will  significantly 
influence  this  important  device  parameter.  Nonetheless, 
these  calculations  do  offer  substantial  insight  into  pro- 
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blems  inherent  to  the  design,  development,  and  fabrica- 
2 

tion  of  I L type  semiconductor  devices. 

4 . 0 Conclusions 

Numerous  conclusions  can  be  drawn  from  the  foregoing 
2 

calculations  of  I L vertical  transistor  operation.  These 
conclusions  can  be  directed  toward  many  aspects  of  modeling 
this  semiconductor  device;  toward  problems  associated  with 

device  fabrication;  and  toward  additional  model  develop- 
2 

ments  for  I L operation.  It  is  suggested  that  these  con- 
clusions offer  a greater  level  of  accuracy  than  would  be 
realized  for  other,  new,  semiconductor  devices;  this  in- 

2 

creased  accuracy  arises  from  the  similarity  between  an  I L 

vertical  transistor  and  the  conventional  bipolar  transis- 

2 

tor.  Many  of  the  calculated  properties  for  I L vertical 
transistor  operation  have  been  previously  encountered 
during  mathematical  studies  of  the  bipolar  transistor. 

From  this  analysis,  in  conjunction  with  similar  studies 
for  bipolar  operation,  it  is  suggested  that  the  I L tech- 
nology could  introduce  a greater  level  of  engineering  dif- 
ficulty than  the  bipolar  technology.  It  is  also  suggested 
that  this  increased  level  of  difficulty  could  be  particu- 
lar y evident  in  the  design  and  development  of  high-speed 

structures,  containing  a narrow  base  width.  During  down 
2 

opei ation  the  I L vertical  transistor  is  similar  to  a con- 
-/en  ional  b polar  transistor.  For  this  reason,  there  is 
xi.t  le  reason  to  expect  a significant  simplification  of 
the  known  problems  associated  with  bipolar  design,  develop- 
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ment,  and  fabrication.  In  addition  to  these  problems, 
many  additional  difficulties  appear  evident  when  con- 

2 

sidering  the  required  operating  characteristics  for  I L 
up  operation. 

The  foregoing  analysis  implies  that  formidable  pro- 
blems could  arise  from  the  small  "built-in"  voltage 
associated  with  an  up  emitter  junction.  This  "built-in" 
voltage  is  substantially  less  than  that  encountered  at 
the  down  emitter  and,  in  addition,  the  up  emitter  area  is 
substantially  greater  than  that  of  the  down  emitter.  In 
combination,  these  inherent  properties  could  produce  im- 
portant limitations  in  the  electric  current  obtainable  from 
a forward  biased  up  emitter.  An  excessive  forward  bias  of 
the  up  emitter,  relative  to  its  "built-in"  voltage,  could 
produce  a reduced  up  emitter  injection  efficiency,  and  an 
excessive  emitter  junction  space-charge  layer  capacitance. 

Conclusions  drawn  from  current  gain  calculations  for 
up  operation  are  consistent  with  other  published  investiga- 
tions of  this  subject,  both  experimental  and  theoretical. 

These  calculations  indicate  the  up  operation  current  gain 
2 

of  an  I L structure  will  be  greatly  influenced  by  the 
distance  between  its  up  emitter  and  its  buried  layer. 

These  calculations  indicate  that  a relatively  small  in- 
crease of  distance  between  the  up  emitter  and  buried  layer 
could  produce  a substantial  decrease  of  current  gain. 

As  in  bipolar  transistor  operation,  the  up  emitter 
to  buried  layer  distance  also  influences  the  reach-through 
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voltage  during  down  operation,  assuming  space-charge  layer 

interference.  Therein  a basic  problem  could  be  encountered 
2 

in  I L design.  The  need  for  a close  proximity  between  the 
up-emitter  and  the  buried  layer  (for  up  emitter  injection 
efficiency)  could  produce  an  excessively  small  down  opera- 
tion reach-through  voltage.  Clearly,  an  extensive  analytical 

2 

study  might  be  required  during  I L design  as  a means  to  ob- 
tain the  necessary  "trade-off"  for  satisfactory  device 
operation. 

Unlike  the  bipolar  transistor,  during  up  operation 

2 

an  I L vertical  transistor  exhibits  large  changes  of  elec- 
trical base  width,  with  changes  of  emitter  junction  biasing 
voltage.  This  situation  is  attributable  to  the  relatively 
small  impurity  atom  density  gradient  associated  with  this 
junction.  A small  p-n  junction  impurity  atom  density  gra- 
dient implies  a relatively  large  equilibrium  space-charge 
layer  width;  hence,  forward  biasing  the  up  emitter  junc- 
tion produces  a reduction  of  space-charge  layer  width,  and 
an  associated  increase  of  electrical  base  width. 

It  is  also  suggested  that  an  excessively  small  open- 

base  breakdown  voltage  (BV,,__)  could  represent  an  inherent 

2 

difficulty  associated  with  I L up  operation.  Although  this 
parameter  has  not,  as  yet,  been  included  in  the  present 
model,  extensive  understanding  of  this  problem  has  been 
gained  from  previous  studies  of  bipolar  transistor  opera- 
tion. Briefly,  open-base  breakdown  arises  when  carrier 
multiplication  in  the  up  collector,  M,  becomes  excessive; 

instability  usually  occurs  if 
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approximately.  In  bipolar  transistor  operation,  Eq.  44  is 
usually  satisfied  when  the  applied  collector  junction 
biasing  voltage  is  about  24%  to  28%  of  the  collector  junc- 
tion avalanche  breakdown  voltage. 

The  up  operation  collector  junction  usually  exhibits  a 
breakdown  voltage  of  about  4.5  volts.  As  a consequence, 
an  up  collector  biasing  voltage  of  2.0  volts  could  produce 
a substantial  amount  of  carrier  multiplication  and,  hence, 

a reduced  BV  . For  example,  Eq.  44  indicates  that  an  M 

2 

of  only  1.5  would  produce  instability  in  an  I L device 

exhibiting  a current  gain  of  2;  this  magnitude  of  M could 

2 

be  very  close  to  the  operating  range  for  an  I L structure. 

2 

This  model  for  I L operation  has  been  used  to  emphasize 
the  calculated  properties  during  up  operation.  This  model 
remains  incomplete,  yet  in  this  incomplete  stage  of  develop- 
ment substantial  insight  can  be  gained  concerning  problems 
2 

of  I L design  and  development.  Because  the  model  discussed 

here  was  drawn  from  previous  studies  of  bipolar  transistor 

operation,  it  is  equally  applicable  to  a similar  study  of 

2 

down  operation  for  an  I L structure.  It  is  suggested  that 
all  work  necessary  for  such  an  extension  could  be  completed 
in  about  two  months  of  effort. 

Including  in  this  model  the  mechanisms  associated  with 
open-base  breakdown  (both  up  and  down  operation)  represents 
a relatively  trivial  task.  The  computational  tools  needed 


to  calculate  space-charge  layer  multiplication  are  already 
operational;  one  need  only  couple  this  calculation  with 
the  current  gain. 

As  previously  stated,  during  up  operation  the  base 
region  potential  distribution  is  known  to  produce  lateral 
variations  of  up  emitter  biasing  voltage.  This  problem 
was  previously  encountered  in  the  analysis  of  bipolar  tran- 
sistor operation:  during  studies  of  emitter  crowding  and 

collector  junction  saturation.  A consequence  of  this 

previous  experience  is  a knowledge  of  how  to  implement  such 

2 

mechanisms  into  a model  of  I L operation.  Briefly,  it  is 
suggested  that  by  coupling  through  base  resistances  one- 
dimensional  lumped  models  for  each  intrinsic  I L region, 
an  adequate  model  could  be  obtained  for  variations  of  cur- 
rent gain  between  each  up  collector  junction. 

Although  complicated  in  appearance,  this  proposed  model 
is  computationally  simple  and  requires  very  little  computer 
time.  In  reality,  the  entire  model  has  been  implemented  on 
a small,  programmable,  desk  calculator.  Through  previous 
experience  with  this  model,  in  conjunction  with  a knowledge 
of  relative  computer  speeds,  it  is  estimated  that  the  present 
model  would  require  about  15-20  seconds  CPU  time  on  a large 
computer  facility  (typically  IBM  360/85) . 

It  is  also  suggested  that  this  model  is  based  upon  im- 
portant physical  mechanisms  associated  with  semiconductor 
device  operation.  For  this  reason,  introducing  into 


this  model  mechanisms  associated  with  radiation  damage  would 
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not  represent  an  unreasonable  engineering  problem.  In 


combination,  this  model  would  offer  an  ability  to  inex- 
pensively calculate  the  electrical  and  radiation  proper- 
ties  of  an  I L vertical  transistor.  It  is  suggested  that 

^ a computational  tool  of  this  type  could  be  completed  and 

y 


made  available  to  the  semiconductor  industry  with  about 
k one  additional  year  of  development  effort. 


Glossary  of  Symbols 


BV 

CEO 

a 

6 

C(x) 

C 


'B 

'1,2 

'eq 


D 


n 


n 


ri 


k 

K 

K 

L 

M 

y 

n 


n 


n . 

1 

P 

q 


ext 

^int 


T 

T 

r 

V 


av 


Open  base  breakdown  voltage 
Common  base  current  gain 
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Impurity  concentration  profile 
Surface  impurity  concentration 
Epitaxial  layer  impurity  concentration 
Junction  capacitances 
Equilibrium  junction  capacitance 
Electron  diffusion  coefficient 

Average  electron  diffusion  coefficient  in  base  region 

Electric  field  intensity 

Permittivity  of  free  space 

Collector  current  density 

Electron  current  density 
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Relative  permittivity  of  silicon 
Impurity  atom  diffusion  length 
Carrier  multiplication  factor 
Electron  mobility 
Electron  density 
Intrinsic  density 
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Electrostatic  potential 
Electron  charge  (magnitude) 

Stored  minority  charge  in  the  "up"  emitter  region 
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Stored  minority  charge  in  the  intrinsic  base  region 
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Electron  lifetime  in  the  base  region 
Applied  reverse  bias  voltage 
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Collector-base  bias  voltage 
Emitter-base  bias  voltage 
Equilibrium  junction  voltage 

Reach  through  voltage,  "up"  emitter  unbiased 

Reach  through  voltage,  "up"  emitter  forward  biased 
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DEPARTMENT  OF  THE  NAVY 
ARLINGTON,  VA  22217 
ATTN  TECHNICAL  LIBRARY 
+ATTN  CODE  1427 

ATTN  CODE  I46U,  THOMAS  P.  QUINN 
ATTN  DR.  DAVID  LEWIS 

COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
HEADQUARTERS 
WASHINGTON,  DC  21360 
ATTN  TECH  LIB 

ATTN  LCDR  HUGO  HARDT  (AIR-350-F) 


COMMANDER 

NAVAL  SEA  SYSTEMS  COMMAND 
NAVY  DEPARIMENT 
WASHINGTON,  DC  20362 
+ATTN  sea-9931,  RILEY  B.  LANE 
*ATTN  sea-9931,  SAMUEL  A.  BARHAM 

COMMANDER 

NAVAL  SHIP  ENGINEERING  CENTER 

CENTER  BUILDING 

HYATTSVILLE,  MD  20782 

+ATTN  CODE  617UD2,  EDWARD  F.  DUFFY 


US  ARMY  ORDNANCE  CENTER  & SCHOOL 
ABERDEEN  PROVING  GROUND,  MD  21005 
ATTN  ATSI-CTD 

US  ARMY  INTELLIGENCE  CENTER  & SCHOOL 
FT.  HUACHUCA,  AZ  856 13 
ATTN  ATSI-CTD 

COMDT 

US  ARMY  ENGINEER  SCHOOL 
FT.  BELVOIR,  VA  22060 
ATTN  ATSE-CDT 


COI^IANDER 

NAVAL  WEAPONS  CENTER 

CHINA  LAKE,  CA  93555 

+ATTN  CODE  533,  TECHNICAL  LIBRARY 

COMMANDING  OFFICER 
NAVAL  WEAPONS  EVALUATION  FACILITY 
KIRTLAND  AIR  FORCE  BASE 
ALBUQUERQUE,  NM  87117 
ATTN  LAWRENCE  R.  OLIVER 
+ATTN  CODE  ATG,  MR.  STANLEY 


US  ARMY  INFANTRY  SCHOOL 
FT.  BENNING,  GA  31905 
ATTN  ATSA-CD 

COMMANDING  OFFICER 
NAVAL  INTELLIGENCE  SUPPORT  CENTER 
I43OI  SUITLAND  ROAD,  BLDG  5 
WASHINGTON,  DC  20390 
ATTN  TECHNICAL  LIBRARY 


COMMANDER 

NAVAL  SURFACE  WEAPONS  CENTER 
DAHLGREN  LABORATORY 
DAHLGREN,  VA  22I4I48 

ATTN  CODE  FUR,  ROBERT  A.  AMADORI 

COMMANDER 


SUPERINTENDENT 

NAVAL  POSTGRADUATE  SCHOOL 

MONTEREY,  CA  939l*0 

ATTN  CODE  212U,  TECH  RPTS  LIBRARIAN 


NAVAL  TELECOMMUNICATIONS  COMMAND 
NAV  TEL  COM  HEADQUARTERS 
I4I4OI  MASSACHUSETTS  AVE,  NW 
WASHINGTON,  DC  20390 
ATTN  TECH  LIB 
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INSTITUTE  FOR  DEFENSE  ANALYSES 
1*00  ARMY-NAVY  DRIVE 
ARLINGTON,  VA  22202 

ATTN  10,  LIBRARIAN,  RUTH  S.  SMITH 

IRT  CORPORATION 
7650  CONVOY  COURT 
SAN  DIEGO,  CA  92138 
ATTN  R.  L.  MERTZ 
ATTN  ERIC  P.  WENAAS 
ATTN  RALPH  H.  STAHL 
ATTN  J.  AZAREWICZ 
ATTN  JAMES  NABOR 
ATTN  T.  A.  TUMOLILLO 

JET  PROPULSION  LAB 
CALIFORNIA  INSTITUTE  OF  TECHNOLOGY 
1*800  OAK  GROVE  DR 
PASADENA,  CA  91103 
ATTN  ALAN  G.  STANLEY 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  208l0 
*ATTN  PETER  E.  PARTRIDGE 


KAMAN  SCIENCES  CORPORATION 
P.O.  BOX  71*63 

COLORADO  SPRINGS,  CO  80933 
+ATTN  W.  FOSTER  RICH 
+ATTN  WALTER  E.  WARE 

LITTON  SYSTEMS,  INC. 

GUIDANCE  & CONTROL  SYSTEMS  DIVISION 
5500  CANOGA  AVENUE 
WOODLAND  HILLS,  CA  91361* 

+ATTN  VAL  J.  ASHBY,  MB  67 

NORTHROP  CORPORATION 
ELECTRONIC  DIVISION 
1 RESEARCH  PARK 

PALOS  VERDES  PENINSULA,  CA  92O7I* 
*ATTN  BOYCE  T.  AHLPORT 


MARTIN  MARIETTA  AEROSPACE 

ORLANDO  DIVISION 

P.O.  BOX  5837 

ORLANDO,  FL  32805 

+ATTN  WILLIAM  W.  MRAS,  MP-1*13 

martin  MARIETTA  CORPORATION 

DENVER  DIVISION 

PO  BOX  179 

DENVER,  CO  80201 

+ATTN  J.  E.  GOODWIN,  MAIL  Ol*52 

MCDONNELL  DOUGLAS  CORPORATION 
POST  OFFICE  BOX  5l6 
ST.  LOUIS,  MO  63166 
+ATTN  TOM  ENDER 

MCDONNELL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  9261*7 
ATTN  W.  R.  SPARK,  13-3 

MCDONNELL  DOUGLAS  CORPORATION 
PO  BOX  1850 

ALBUQUERQUE,  NM  87103 
+ATTN  THOMAS  J.  LUNDREGAN 

GEORGE  C.  MESSENGER 
1501  SERENADE  TERRACE 
CORONA  del  MAR,  CA  92625 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
ATTN  CONRAD  L.  LONGMIRE 

mitre  corporation,  the 

ROUTE  62  AND  MIDDLESEX  TURNPIKE 
P.O.  BOX  200 
BEDFORD,  MA  01730 
+ATTN  M.  F.  FITZGERALD 

R & D ASSOCIATES 
PO  BOX  9695 

MARINA  DEL  REY,  CA  90291 
+ATTN  S.  CLAY  ROGERS 


! ( 
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COMMMDJi'jK 

USA  ARMAMENTS  COMMAND 
ROCK  ISLAND,  IL  61201 
ATTN  DRSAR-ASF,  FUZE  DIV 
ATTN  DRSAR-RDF,  SYS  DEV  DIV  - FUZES 
ATTN  DRSAR-RDS 


DEFENSE  DOCUMENTATION  CENTER 
CAMERON  STATION  BLDG  5 
ALEXANDRIA,  VA  223ll4 
+ATTN  TC-TCA,  (12  COPIES) 


OFC,  CHIEF  OF  RESEARCH  & DEVELOPMENT 
USA  RESEARCH  & DEV  GROUP  (EUROPE) 

BOX  15 

FPR  NEW  YORK  09510 
ATTN  LTC  EDWARD  E,  CHICK 

CHIEF,  MATERIALS  BRANCH 


NATIONAL  COMMUNICATIONS  SYSTEM 
OFFICE  OF  THE  MANAGER 
WASHINGTON,  DC  20305 

ATTN  NCS-TS,  CHARLES  D.  BODSON 


DIRECTOR 

NATIONAL  SECURITY  AGENCY 

COMMANDER  FT.  GEORGE  G.  MEADE,  MD  20755 

USA  MISSILE  & MUNITIONS  CENTER  & SCHOOL  +ATTN  0,  0.  VAN  C-UNTEN-R-U25 
REDSTONE  ARSENAL,  AL  35809  +ATTN  D.  VINCENT 

ATTN  ATSK-CTD-F  +ATTN  P.  DeBOY 

ATTN  ATSK-CTD  +ATTN  E,  BUTALA 

R233/LTC  J.  BOHREN 


director 

DEFENSE  INTELLIGENCE  AGENCY 
WASHINGTON,  DC  20301 
+ATTN  DI-TD,  MR.  EDWARD  OFARRELL 
ATTN  TECHNICAL  LIBRARY 


DIRECTOR 

DEFENSE  COMMUNICATIONS  AGENCY 
WASHINGTON,  DC  20305 
ATTN  CODE  540,5 
ATTN  CODE  1+30 

+ATTN  CODE  930,  MONTE  I.  BURGETT,  FR 


DEFENSE  COMMUNICATION  ENGINEER  CENTER 
i860  WIEHLE  AVENUE 
RESTON,  VA  22090 
ATTN  CODE  R720,  C.  STANSBERRY 
+ATTN  CODE  Rl+10,  JAMES  W.  MCLEAN 
»ATTN  CODE  R320,  C.  W.  BERGMAN 
ATTN  CODE  Rl+00 


DIRECTOR 

INTERSERVICE  NUCLEAR  WEAPONS  SCHOOL 
KIRTLAND  AFB,  NM  87115 
ATTN  TECH  LIB 

COMMANDER 
FIELD  COMMAND 
DEFENSE  NUCLEAR  AGENCY 
KIRTLAND  AFB,  NM  87115 
+ATTN  FCPR 

ATTN  FCSM-F3  CDR  SMITH 

DIR  OF  DEFENSE  RSCH  & ENGINEERING 
DEPARTMENT  OF  DEFENSE 
WASHINGTON,  DC  20301 
ATTN  DD/S&SS 


OJCS /J-6 
THE  PENTAGON 
WASHINGTON,  DC  20301 
ATTN  J-6,  ESD-2 

*ATTN  J-3  RDTA  BR  WWMCCS  PLANS  DIV 
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COMMANDER  IN  CHIEF 
STRATEGIC  AIR  COMMAND 
OFFUTT  AFB,  NB  68113 
+ATTN  NRI-STINFO  LIBRARY 

DIVISION  OF  MILITARY  APPLICATION 
US  ENERGY  RSCH  & DEV  ADMIN 
WASHINGTON,  DC  205^*5 
ATTN  TECHILTCAL  LIBRARY 

EG&G,  INCc 

LOS  ALAMOS  DIVISION 

PO  BOX  809 

LOS  ALAMOS,  NM  85514*4 
ATTN  L.  DETCH 

LOS  ALAMOS  SCIENTIFIC  LABORATORY 
P.O.  BOX  1663 
LOS  ALAMOS,  NM  875*414 
ATTN  ARTHUR.  FREED 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
PO  BOX  969 

LIVERMORE,  CA  9*4550 
ATTN  TECHNICAL  LIBRARY 

SANDIA  LABORATORIES 
PO  BOX  5800 

ALBUQUERQUE,  NM  8TH5 
ATTN  ORD  9353,  R.  L.  PARKER 

US  ENERGY  RSCH  & DEV  ADMIN 
ALBUQUERQUE  OPERATIONS  OFFICE 
PO  BOX  5*400 

ALBUQUERQUE,  NM  87115 
+ATTN  WSSB 


UNION  CAP.BIDE  CORPORATION 
HOLIFIELD  NATIONAL  LABORATORY 
P.O.  BOX  X 

OAK  RIDGE,  TN  37830 
ATTN  PAUL  R.  BARNES 


SAMSO/DY 

POST  OFFICE  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
+ATTN  DYS,  MAJ  LARRY  A,  DARDA 

SAMSO/MN 

NORTON  AFB,  CA  92*.09 
ATTN  MNNIi,  CAPT  B.  STEWART 
ATTN  MNini,  CAPT  MICHAEL  V.  BELL 
*ATTN  MNNG,  CAPT  DAVID  J.  STRCBFI 

SAJ.ISO/SK 

POST  OFFICE  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
+ATTN  SKF,  PETER  K.  STADLER 

SAMS 0 /YD 

POST  OFFICE  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CA  90009 
+ATTN  YDD,  MAJ  M.  F.  SCHNEIDER 

CORIANDER  IN  CHIEF 
STRATEGIC  AIR  COMMAND 
OFFUTT  AFB,  NB  68113 
+ATTN  NRI-STINFO  LIBRARY 

DIVISION  OF  MILITARY  APPLICATION 
US  ENERGY  RSCH  & DEV  AIMIN 
WASHINGTON,  DC  205*45 
ATTN  TECHNICAL  LIBRARY 

EG&G,  INC. 

LOS  ALAMOS  DIVISION 
PO  BOX  809 

LOS  ALAMOS,  NM  855*4*4 
ATTN  L.  DETCH 

LOS  ALAMOS  SCIENTIFIC  LABORATORY 
P.O.  BOX  1663 
LOS  ALAMOS,  NM  875*4*4 
ATTN  ARTHUR  FREED 
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UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
PO  BOX  808 
LIVERMORE,  CA  91*550 
+ATTN  HANS  KRUGER,  L-96 


SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
PO  BOX  969 

LIVERMORE,  CA  91*550 
ATTN  TECHNICAL  LIBRARY 


CENTRAL  INTELLIGENCE  AGENCY 
ATTN:  RD/SI  RM  5G1+8,  HQ  BLDG 
WASHINGTON,  DC  20505 

HQ  USAF/RD 

WASHINGTON,  DC  20330 
ATTN  RDQPN 


HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION,  (AFSC) 
L,  G.  HANS COM  FIELD 
BEDFORD,  MA  01730 
+ATTN  YWEI 

+ATTN  YSEV,  LTC  DAVID  C.  SPARKS 
»ATTN  DCKE  LEWIS  STAPLES 
ATTN  TECHNICAL  LIBRARY 
ATTN  XRT,  LTC  JOHN  M„  JASINSKI 


COMMANDER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT-PATTERSON  AFB,  OH  1*51+33 
ATTN  TD-BTA,  LIBRARY 


CENTRAL  INTELLIGENCE  AGENCY 
ATTN:  RD/SI  RM  5Gl*8,  HQ  BLDG 
WASHINGTON,  DC  2o505 
ATTN  Tl'.CHNICAL  LIBRARY 


SANDIA  LABORATORIES 
PO  BOX  5800 

ALBUQUERQUE,  NM  87115 
ATTN  ORD  9353,  R.  L.  PARKER 


US  ENERGY  RSCH  & DEV  ADMIN 
ALBUQUERQUE  OPERATIONS  OFFICE 
PO  BOX  51*00 

ALBUQUERQUE,  NM  87115 
+ATTN  WSSB 

UNION  CARBIDE  CORPORATION 
HOLIFIELD  NATIONAL  LABORATORY 
P.O.  BOX  X 
OAK  RIDGE,  TN  37830 
ATTN  PAUL  R.  BARNES 


UNIVERSITY  OF  CALIFORNIA 
LAV7RENCE  LIVERMORE  LABORATORY 
PO  BOX  808 
LIVERMORE,  CA  91*550 
+ATTN  HANS  KRUGER,  L-96 


AEROSPACE  CORPORATION 
PO  BOX  92957 
LOS  ANGELES,  CA  9OOO9 
+ATTN  JULIAN  REINHEIMER 


AEROJET  ELECTRO-SYSTEMS  CO.  DIV. 
AEROJET-GENERAL  CORPORATION 
P.O.  BOX  296 
AZUSA,  CA  91702 

+flTTN  THOMAS  D.  HANSCOME.  B170/D6711 


BATTELLE  IffiMORIAL  INSTITUTE 
505  KING  AVENUE 
COLUI>IBUS,  OH  1*3201 
ATTN  STOIAC 
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AERONUTRONIC  FORD  CORPORATION 
AEROSPACE  & COMMUNICATIONS  OPS 
AERONUTRONIC  DIVISION 
FORD  & JAMBOREE  ROADS 
NIWPORT  BEACH,  CA  92663 
+ATTN  KEN  C.  ATTINGER 


HONEYWELL  INCORPORATED 
GOVERNMENT  AND  AERONAUTICAI. 
PRODUCTS  DIVISION 
2600  RIDGEWAY  PARKWAY 
MINNEAPOLIS,  M 55^13 
+ATTN  RONALD  R.  JOHNSON,  Al622 


BENDIX  corporation,  THE 

guidance  systems  division 

TETERBORO,  NJ  OT6o8 
ATTN  TECH  LIB 


GENERAL  ELECTRIC  COMPANY 
TET'IPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  ( PO  DRAWER  QQ) 
SANTA  BARBARA,  CA  93102 
+ATTN  DASIAC,  M.  ESPIG 


BOEING  COMPANY,  THE 

PO  BOX  3707 

SEATTLE,  WA  9 812 >4 

+ATTN  HOWARD  W.  WICKLEIN,  IE  17-H 

»ATTN  ROBERT  S„  CALDWELL.  2R-00 


HONEYWELL  INCORPORATED 

AEROSPACE  DIVISION 

13350  US  HIGHWAY  19 

ST.  PETERSBURG,  FL  33733 

+ATTN  HARRISON  H.  NOBLE,  IE  725-5A 


BROWN  ENGINEERING  COMPANY,  INC. 
CUMMINGS  RESEARCH  PARK 
HUNTSVILLE,  AL  35807 
ATTN  TECH  LIB,  MS12,  P.  SHELTON 

BURROUGHS  CORPORATION 
FEDERAL  AND  SPECIAL  SYSTENB  GROUP 
CENTRAL  AVE  AND  ROUTE  252 
PO  BOX  517 
PAOLI,  PA  29301 
ATTN  ANGELO  J.  MAURIELLO 


GENERAL  ELECTRIC  COMPANY 
RE-ENTRY  S:  ENVIRONMENTAL  SYSTE14S  DIV 
PO  BOX  7722 
3198  CHESTNUT  STREET 
PHILADELPHIA,  PA  I9IOI 
+ATTN  JOHN  W.  PALCHEFSHY,  JR. 

HUGHES  AIRCRAFT  COMPANY 
CENTINELLA  AVENUE  & TEALE  STREETS 
CULVER  CITY,  CA  90230 
+ATTN  M.S.  DI57,  KEN  WALKER 


CHARLES  STARK  DRAPER  LABORATORY  INC. 
68  ALBAirY  STREET 
CAMBRIDGE,  M 02139 
+ATTN  KEiraETH  FERTIG 


GRUMMAN  AEROSPACE  CORPORATION 
SOUTH  OYSTER  BAY  ROAD 
BETHPAGE,  NY  II71I 
+ATTN  JERRY  ROGERS,  DEPT  533 


GENERAL  ELECTRIC  COI^IPAITY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
P.O.  BOX  8555 
PHILADELPHIA,  PA  I9IOI 
+ATTN  JAMES  P.  SPRATT 
ATTN  DANIEL  EDELMAN 
ATTN  DANTE  M.  TASCA 


HUGHES  AIRCRAFT  COMPANY 

SPACE  SYSTEMS  DIVISION 

P.O.  BOX  92919 

LOS  ANGELES,  CA  90009 

+ATTN  WILLIAM  W.  SCOTT,  ME  AIO8O 
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PROJECT  MANAGER 
US  ARMY  TACTICAL  DATA  SYSTEMS 
FORT  MONMOUTH,  NJ  07703 
ATTN  TECH  LIBRARY 
+ATTN  DWAINE  B.  HUEWE 
ATTN  *DRCPN-TDS-SD 
ATTN  DRCPN-TDS/MR.  S.  PSLOSI 


PM 

SANG 

5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22333 
ATTN  DRCPM-NG 

RM 

DRAGON 

US  ARMY  MISSILE  COMMAND 
REDSTONE  ARSENAL,  AL  35809 

PM 

HELLFIRE  MISSILE  SYSTEM 
US  ARMY  MISSILE  COMMAND 
REDSTONE  ARSENAL,  AL  35809 

PM 

HIGH  ENERGY  LA'  R SYSTEM 
US  ARMY  MISSILE  COMMAND 
REDSTONE  ARSENAL,  AL  35809 

PM 

KUWAIT 

US  ARm  MISSILE  COMMAND 
REDSTONE  ARSENAL,  AL  35809 

IW 

CANNON  ARTILLERY  WEAPONS  SYSTEMS 
US  ARMY  ARMAMENT  COMMAND 
ROCK  ISLAND,  IL  61201 

PROGRAM  MANAGER 

Xll-l  TANK  SYS 

28150  DeQUIDRE 

NATICK,  MI  01760 

ATTN  DRCPM-GCM-SW/f® . WOOLCOT 

ATTN  DRCPM-GCM-SW/R.  SLAUGHTER 


PM 

SAFEGUARD  MUNITIONS 
PICATINNY  ARSENAL 
DOVER,  NJ  07801 

PM 

SELECTED  AMMUIMTION 
PICATINNY  ARSENAL 
DOVER,  NJ  07801 

COMMANDER 

US  ARMY  TEST  AND  EVALUATION  COMMAND 
ABERDEEN  PROVING  GROUND,  MD  21005 
+ATTN  AMSTE-EL,  R.  I.  KOLCHIN 
+ATTN  AMSTE-NB,  R.  R.  GALASSO 
ATTN  TECHNICAL  LIBRARY 
ATTN  COL  GREGORY 


COMMANDER 

WHITE  SANDS  MISSILE  RANGE 
WHITE  SANDS  MISSILE  RANGE,  NM  88002 
ATTN  TECHNICAL  LIBRARY 
+ATTN  STEWS-TE-NT 
ATTN  STEWS-TE-AN/MR.  A DeLaPAZ 

COMMANDER 

US  ARMY  SECURITY  AGENCY 
ARLINGTON  HALL  STATION 
UOOO  ARLINGTON  BLVD 
ARLINGTON,  VA  22212 
+ATTN  lARD-T,  DR.  R.  H.  BURKHARDT 
ATTN  TECHILECAL  LIBRARY 

PM 

ADVANCED  ATTACK  HELICOPTER 
12TH  & SPRUCE  STS 
ST  LOUIS,  MO  63160 

ATTN  DRCPM-AAH/MR  ROMANO 

US  AMC  INFANTRY  R&D  LIAISON  OFFICE 
DEPARTMENT  OF  THE  ARMY 
FT.  EENNING,  GA  31905 
ATTN  DRCLR  ( COL  TULLY) 

CRD 

US  ARMY  CONCEPTS  ANALYSIS  AGENCY 
BETHESDA,  MD  200114 
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COMMANDER 

US  ARMY  NATICK  LAB 
NATICK,  MA  01760 
ATTN  STSNL-WR/MR  T.  MARTIN 

COMMANDER 

US  AROT  NUCLEAR  AGENCY 
FORT  BLISS,  TX  79916 
ATTN  TECH  LIB 

ATTN  ATCN-W  COL  DEVERILL  (2  COPIES) 
COMMANDER 

ROCK  ISLAND  ARSENAL 
ROCK  ISLAND,  IL  61201 
ATTN  SARRI-LR/MR  W.  McGARVEY 

COMMANDER 

US  ARMY  SATELLITE  COMM  AGENCY 
FT.  MONMOUTH,  NJ  07'^03 
ATTN  DRCPM-5C-SB/MR  P.  KENNEDY 
ATTN  DRCPM-SA/MR.  MARES CA 
ATTN  DRCPM-SA/V.  CHIWEY 

PROJ  MGR 
PERSHING  SYSTEM 
REDSTONE  ARSENAL,  AL  35809 
ATTN  DRCPM- PE/MR.  J.  PETTIT 

COMMANDEii 

US  ARMY  MISSILE  COMMAND 
REDSTONE  ARSENAL 
HUNTSVILLE,  AL 
ATTN  DPSMI-RGP, 

+ATTN  DRCPM- LC EX, 

ATTN  DRSMI-RRR/DR. 

COMMANDER 
US  ARMY  TANK  AUTOMOTIVE  C0M4AND 
WARREN,  MI  U8O89 

+ATTN  AMSTA-RHM,  ILT  PETER  A.  HASEK 
ATTN  TECH  LIBRARY 

-'MMANDER 

: APf'*Y  MATERIALS  & MECHANICS 
' r-rFAPCH  CENTER 
"-F-T^wri,  MA  02172 

.-  '"S  'PXMP-HH,  JOHN  F.  DIGNAM 


COMMANDANT 

US  ARMY  FIELD  ARTILLERY  SCHOOL 
FORT  SILL,  OK  73503 
*ATTI^  ATSFA-CTD-ME,  HARLEY  MOBERG 
ATTN  TECH  LIBRARY 
ATTN  ATSF-CD 

US  ARMY  INFANTRY  SCHOOL 
NBC  COMMITTEE,  BROD 
FT.  BENNING,  GA  31905 
ATTN  ASH-B 

PROJ  MGR,  SAM-D 
REDSTONE  ARSENAL,  AL  35809 
ATTN  DRCPM-IffiDI/I'«.  R.  HASE 

PROJ  MGR 

US  ARMY  TACTICAL  COMM  SYS 
FT.  MONMOUTH,  NJ  07703 
ATTN  DRCPM-ATC/COL  DOBBINS 

PROJ  MGR 

AVIATION  SYSTEMS  INTEGRATION 
USA  AVIATION  SYSTEMS  COMMAND 
STo  LOUIS,  MO  63160 

PROJ  MGR 

NAVIGATION  CONTROL  SYSTE^e 
USA  ELECTRONICS  COMMAND 
BLDG  2525 

FT.  MONMOUTH,  NJ  07703 
CRD 

USA  PROCUREMENT  RESEARCH  OFFICE 
LEE,  VA  23801 

DEFENSE  SYSTEMS  MJU^AGEMEOT  SCHOOL 
FT.  BELVOIR,  VA  22060 
ATTN  MR.  PURNELL 

PM 

FAMECE 

FT.  BELVOIR,  VA  22060 
PROJ  MGR 

MULTI -SERVICE  COMMUNICATIONS  SYSTEMS 
USA  ELECTRONICS  COIWlND 
FT.  MONMOUTH,  NJ  07703 


35809 

HUGH  GREEN 
HOWARD  H.  HENRIKSEIT 
F.  GIBSON 


NAVAL  ELECTRONIC  SYSTEMS  COMMAND 
HEADQUARTERS 
WASHINGTON,  DC  20360 
AT'-N  TECH  LIB 

ATTN  PME11T-215A,  GUNTER  BRUNHART 
+ATTN  PME  117-21 
ATTN  PME  IIT-T 

*ATTN  CODE  5032  CHAS . W.  NEILL 
ATTN  CODE  503/C  FRANK  BERG 

DIRECTOR 

NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  DC  20 375 

ATTN  CODE  6631,  JAMES  C.  RITTER 
ATTN  CODE  UOOU,  EMANUAL  L.  BRAI^CATO 
ATTN  CODE  2027,  TECH  LIB 
ATTN  CODE  7701,  JACK  D.  BROWN 
ATTN  CODE  7706,  JAY  P.  BORIS 
ATTN  CODE  h6h,  R,  GRACEN  JOINER 
ATTN  CODE  6600 , B.  FARADAY 
ATTN  CODE  6600 , E.  WOLICKI 
' ATTN  CODE  6600 , N,  WILSEY 
ATTN  CODE  6600 , V.  CAROSELLA 
ATTN  CODE  6600 , R.  WENZEL 

COMMANDER 

NAVAL  SURFACE  TOAPONS  CENTER 
VHIITE  OAK,  SILVER  SPRING,  MD  20910 
+ATTN  CODE  730,  TECH  LIB 
+ATTN  CODE  U3I,  EDWIN  B.  DEAN 
+ATTN  CODE  122i*,  NAVY  NUC  PRGMS  OFF 
ATTN  CODE  U3I,  EDWIN  R.  RATHBURN 
+ATTN  CODE  U3I,  JOHN  H.  MALLOY 
ATTN  R.  JENKINS 
ATTN  R.  PARTAK 
ATTN  J,  FRANKLIN 
ATTN  D.  GILCREASE 


COMMANDER 

ADC/DF 

ENT  AFB,  CO  80912 

ATTN  DEEDS,  JOSEPH  C.  BRANNAN 

COMMANDER 

ADC/XP 

ENT  AFB,  CO  80912 

ATTN  XPQD,  MAJ  G.  KUCH 

COMMANDER 

AERONAUTICAL  SYSTEMS  DIVISION,  AFSC 
WRIGHT-PATTERSON  AFB,  OH  U5I+33 
ATTN  TECHNICAL  LIBRARY 
ATTN  ASD-YH-EX 

AF  ARMAMENT  LABORATORY,  AFSC 
EGLIN  AFB,  FL  325i*2 
+ATTN  DLOSL-LIB 

AF  CAMBRIDGE  RSCH  LABS,  AFSC 
L,  Go  HANS COM  FIELD 
BEDFORD,  MA  01730 
ATTN  EDWARD  A,  BURKE 

AF  WEAPONS  LABORATORY,  AFSC 
KIRTLAND  AFB,  NM  87II7 
ATTN  P.  VAIL 
ATTN  J.  NICHOLS 
AFTAC 

PATRICK  AFB,  FL  32925 
ATTN  TECH  LIB 

AIR  FORCE  AVIONICS  LABORATORY,  AFSC 
WRIGHT-PATTERSON  AFB,  OH  i*5U33 

*ATTN  AFAL-TEA  F_ANS  J.  HENNECKE 


DIRECTOR 

STRATEGIC  SYSTEMS  PROJECT  OFFICE 
NAVY  DEPARTMENT 
WASHINGTON,  DC  20376 
+ATTI'f  SP2701,  JOHN  W.  PITSENBERGER 
+ATTN  NSP-23i*2,  RICHARD  L.  COLEMAN 
ATTN  NSP-1<3,  TECH  LIB 
+ATTN  NSP-230,  DAVID  GOLD 
^‘ATTN  NSP-27331,  PHIL  SPECTOR 


HEADQUARTERS 

AIR  FORCE  SYSTEMS  COMMAND 
ANDREWS  AFB 
WASHINGTON,  DC  20331 
ATTN  TECHNICAL  LIBRARY 

COMMANDER-IN-CHIEF 
US  PACIFIC  FLEET 
FPO  SAN  FRANCISCO  966IO 
ATTN  DOCUMENT  CONTROL 
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PROJ  MGR 

REMOTELY  MONITORED  BATTLEFIELD  SENSOR  SYSTEMS 
USA  ELECTRONICS  COMMAND 
FT.  MONMOUTH,  NJ  07703 


CRD 

USA  FOREIGN  SCIENCE  & TECHNOLOGY  CENTER, 
FEDERAL  OFFICE  BLDG 
CHARLOTTESVILLE,  VA  20901 


PM 

SINGLE  CHANNEL  GROUND  AND  AIRBORNE  RADIO  SUBSYSTEM 
US  ARMY  ELECTRONICS  COMMAND 
FT.  MONMOUTH,  NJ  07703 

COMMANDER 

US  ARMY  MOBILITY  EQUIPMENT  R & D CENTER 
FORT  BELVOIR,  VA  22060 

ATTN  STSFB-MW,  JOHN  W.  BOND,  JR. 

ATTN  TECHNICAL  LIBRARY 


COMMANDER 

US  ARMY  TRAINING  AND  DOCTRINE  COMMAND 
FORT  MONROE,  VA  2365I 
ATTN  TECH  LIBRARY 
»ATTN  ATORI-OP-SD 
ATTN  ATCD-PM-S 


COMMANDING  OFFICER 

NUCLEAR  WEAPONS  TRAINING  CENTER  PACIFIC 
NAVAL  AIR  STATION,  NORTH  ISLAND 
SAN  DIEGO,  CA  92135 
ATTN  CODE  50 


WESTINGHOUSE  ELECTRIC  CORPo'rATION 
DEFENSE  AND  ELECTRONIC  SYSTEMS  CENTER 
P.O.  BOX  1693 

FRIENDSHIP  INTERNATIONAL  AIRPORT 
BALTIMORE,  MD  21203 
+ATTN  HENRY  P.  KALAPACA,  MS  3525 


I 


I i 
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COMMMDER 

US  ARMY  MATERIEL  DEVELOPMENT  AND  READINESS  COMMAND 
5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22333 
ATTN  DRCRD,  RES,  DEV,  & ENGR  DIRECTORATE 
ATTN  DRCRD-WB,  JOHN  F.  CORRIGAN 
ATTN  TECHNICAL  LIBRARY 
ATTN  DRCPM 


DIRECTOR 

DEFENSE  NUCLEAR  AGENCY 
WASHINGTON,  DC  20305 

ATTN  PETER  HAAS,  DEP.  DIR,  SCIENTIFIC  TECHNOLOGY 
ATTN  RAEV/CPT  R,  STEAD /CPT  R.  DAVIS  (2  COPIES) 


COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFFISS  AFB,  NY  13l*Ul 
ATTN  EMTLD,  DDC  LIBRARY 
*ATTN  RBAC,  I.  L.  KRULAC 
ATTN  RBRM,  JOSEPH  B.  BRAUER  (l  COPY) 25  COPIES) 
ATTN  RBRM,  CYLDE  LANE 
ATTN  RADC- DDC /EARL  SELOVER 

DIRECTOR 

JOINT  STRATEGIC  TARGET  PLANNING  STAFF,  JCS 
OFFUTT  AFB 
OMAHA,  NB  68113 
*ATTN  JLTW-2 
ATTN  JPST 

ATTN  JSAS,  MAJ  MICHAEL  F.  BARAN 


DEPT  OF  ARMY 

DEPUTY  CHIEF  OF  STAFF  FOR  OPERATIONS  & PLANS 
WASHINGTON,  DC  20310 
ATTN  DAMO-SSN,  NUCLEAR  PLANS  AND  POLICY 
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COMI^AinDER 

NAVAL  UNDERSEA  CENTER 
SAN  DIEGO,  CA  92152 

ATTN  CODE  608,  CLARENCE  F.  RANSTEDT 


COMMANDER 

US  NAVAL  COASTAL  SYSTEMS  LABORATORY 
PANAMA  CITY,  EL  321*01 
ATTN  TECH  LIB 


DIRECTOR 

US  ARMY  BALLISTIC  RESEARCH  LABORATORIES 
ABERDEEN  PROVING  GROUND,  MD  21005 
ATTN  »DRXBR-VL  ROBERT  L.  HARRISON 
ATTN*DRXBR-X  JULIUS  J.  MESZAROS 


COMMANDING  OFFICER 
NAVAL  AMMUNITION  DEPOT 
CRANE,  IN  1*7522 

+ATTN  CODE  7021*,  JAMES  L.  RAMSEY  (l  COPY) 


CHIEF  OF  RES,  DEV  & ACQUISITION 
DEPARTMENT  OF  THE  ARMY 
WASHINGTON,  DC  20310 

ATTN  DAMA-CSM-N,  LTC  G.  T.  OGDEN 


COMMANDER 

NAVAL  ELECTRONICS  LABORATORY  CENTER 
SAN  DIEGO,  CA  92152 

ATTN  CODE  2 1*00,  S.  W.  LICHTMAN 
ATTN  CODE  2200  1,  VERNE  E.  HILDEBRAND 
ATTN  CODE  3100,  E.  E.  MCCOWN 

US  ARMY  INSTITUTE  FOR  MILITARY  ASSISTANCE 
FT.  BRAGG,  NC  28307 
ATTN  ATSU-CTD 

Dr.  Peter  Hudson 

Attn;  DRSEL-TL-IS 

U.S.  Army  Electronics  Command 

Fort  Monmouth,  NJ  07703 
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COMMANDER 

ROlffi  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFFISS  AFB,  NY  134Ul 
ATTN  EMTLD,  DDC  LIBRARY 
*ATTN  RBAC,  I.  L.  KRULAC 
ATTN  RBRM,  JOSEPH  B.  BRAUER  (l  COPY) 25  COPIES) 
ATTN  RBRM,  CYLDE  LANE 
ATTN  RADC-DDC/EARL  SELOVER 


FAIRCHILD  CAMERA  AIID  INSTRUMENT  CORPORATION 
U6h  ELLIS  STREET 
MOUNTAIN  VIEW,  CA  9U0UO 

ATTN  2-233,  MR.  DAVID  K.  MYERS 


BETJDIX  CORPORATION,  THE 

RESEARCH  LABORATORIES  DIV 

BENDIX  CENTER 

SOUTHFIELD , MI  U8075 

+ATTN  MGR  PROGM  DEV,  DONALD  J.  NIEHAUS 


AVCO 

GOVERNMENT  PRODUCTS  GROUP 
201  LOWELL  STREET 
WILMINGTON,  MA  OI887 

ATTN  RESEARCH  LIBRARY,  A83O,  RM  7201 


LOCKHEED  MISSILES  AND  SPACE  COMPANY,  INC. 

P.O.  BOX  50»+ 

SUNNYVALE,  CA  91*088 
+ATTN  G.  F.  HEATH,  D/8I-II* 

NATIONAL  ACADEMY  OF  SCIENCES 
2101  CONSTITUTION  AVE,  NW 
WASHINGTON, DC  20Ul8 

+ATTN  DR.  R.  S.  SHANE,  NAT  MATERIALS  ADVISORY  BO 

RAYTHEON  COMPANY 
HARTWELL  ROAD 
BEDFORD,  MA  01730 

+ATTN  GAJANAN  H.  JOSHI,  RADAR  SYS  LAB 


